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ABSTRACT 

H2 pure-rotational emission lines are detected from warm (100-1500 K) molecular gas in 17/55 (31% 
of) radio galaxies at redshift z < 0.22 observed with the Spitzer IR Spectrograph. The summed H2 
0-0 S(0)-S(3) line luminosities are L(H2) = 7 x 10 38 — 2 x 10 42 erg s _1 , yielding warm H2 masses up to 
2 x 1O 1O M0. These radio galaxies, of both FR radio morphological types, help to firmly establish the 
new class of radio-selected molecular hydrogen emission galaxies (radio MOHEGs). MOHEGs have 
extremely large H 2 to 7.7 (im PAH emission ratios: i(H 2 )/L(PAH7.7) = 0.04 — 4, up to a factor 300 
greater than the median value for normal star-forming galaxies. In spite of large H2 masses, MOHEGs 
appear to be inefficient at forming stars, perhaps because the molecular gas is kinematically unsettled 
and turbulent. Low- luminosity mid-IR continuum emission together with low- ionization emission 
line spectra indicate low-luminosity AGNs in all but 3 radio MOHEGs. The AGN X-ray emission 
measured with Chandra is not luminous enough to power the H2 emission from MOHEGs. Nearly 
all radio MOHEGs belong to clusters or close pairs, including 4 cool core clusters (Perseus, Hydra, 
A 2052, and A 2199). We suggest that the H2 in radio MOHEGs is delivered in galaxy collisions or 
cooling flows, then heated by radio jet feedback in the form of kinetic energy dissipation by shocks or 
cosmic rays. 

Subject headings: galaxies: active, galaxies: quasars, galaxies: jets, infrared: galaxies 
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1. INTRODUCTION 

1.1. AGN Radio Jet Feedback 

AGN feedback on the interstellar and intracluster me- 
dia (ISM and ICM) is thought to be an important factor 
regulating gas accretion, star formation, and AGN activ- 
ity in g alaxies dGranatq et al. Ij2~0 04: Sca nnapieco &: OlTI 



20041: IHopkins et al. . 
Ciotti fc Ostriker I [200 



1200 



Croton et al. I [2006; 
Merloni fc Heinz I 120081 ) . 
This may occur through either a sub-relativistic wind 
in radio-quiet quasars or a relativistic radio jet in the 
case of radio-loud quasars and radio galaxies. In the 
latter case, the kinetic energy of the jet integrated over 
its lifetime is typically more than enough to unbind 
the ISM of even the most massive galaxy. While much 
effort has been devote d to simulating the impact of jets 
on the hot ICM fe.g.. lZanni et al.~ll2005l ). less is known 
about the jet-ISM interaction. Recent simulations 
have shown that a radio jet may couple strongly to an 
inhomogeneous, clumpy ISM unt il the radio lobes break 
free into intergalactic space ([Sutherland fc Bicknell I 
2007). In this process, the jet drives radiative shocks 
into dense clumps but does not ablate them, while less 
dense gas may be ablated and advected to larger radii. 
While such detailed simulations have been made of the 
impact of jets on the atomic and ionized components of 
the ISM, the impact on the molecular gas phase has not 
been considered in any detail. 
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lOgle et al. I (|2007f) found that the mid-IR spectrum of 
radio galaxy 3C 326 N displays high-luminosity H2 emis- 
sion, coupled with weak PAH emission indicative of a low 
star formation rate. We presented evidence that 3C 326 
N belongs to an interacting galaxy pair and suggested 
that the copious warm H2 is heated by accretion from 
the companion 3C 326 S. Further investigation revealed 
the presence of fa st neutral gas outflow, possibly driven 
by the radio jet (jNesvadba et al. I [2010V Therefore, a 
strong interaction between the radio jet and molecular 
gas, rather than accretion, may power the H 2 emission. 

H2 emission may prove to be a key in understanding 
how radio jet feedback impacts the molecular gas phase 
and thereby suppresses star formation in massive galax- 
ies. The most likely H2 heating mechanism is jet-driven 
shocks in a multiphase interstellar medium. Cosmic ray 
heating and X-ray emission from the AGN are alternative 
heating mechanisms that we will consider. 

1.2. H2 Emission from Galaxies 

Rotational and ro-vibrational H2 transitions are im- 
portant cooling channels for warm molecular gas, ob- 
served as a series of infrared emission lines. These 
transitions are highly forbidden because of the zero 
permanent dipole moment of H2. In normal star- 
forming galaxies with total H2 to total (3-1100 jum) 
IR ratios L(H 2 )/Ltir < 0.001, these lines are 
thought to be emitted by ste ll ar photodissociation 
regions (iRigopoulou et~al~l l2002t iHigdon et al~l 120061: 
iRoussel et al. II2007D . 

Much stronger molecular hydrogen emission galax- 
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ies (MOHEGs) with L{R 2 )/L Tm > 0.001 - 0.1 have 
been revealed by recent Spitzer observations (jQgle et al. I 
I2007T) . Extreme H 2 emission has been found in a va- 
riety of galaxy types so fa r , inc luding luminous in- 
frared galaxies (Lutz et aT~l 120031 ). colliding galaxies 
(lAppleton et al. I 20061)7 central dominant cluster galax- 



ies (|Egami et al. 1 120061: iDe Messieres et al. II2009I). an d 
radio galaxies ([Ogle et al. 1 120071: iWillett et al. M2010D . 
Enhanced H2 emission has been observed in AGNs 



with ISO (Ri gopoulou et aL~l l2002j ) . and more recently 
from AGNs in the Spitzer SINGS survey (jRoussel et al. I 
|2007[ ). Relatively strong H2 emissio n is also seen in 
PG quasars with strong PAH emission ((Schweitzer et al~l 
[2001 IVeilleux et aD I2009T) . AGN X-ray heating or 
shocks have been cited as possible ways to power the 
H2 emission. 

The first galaxy discovered with strong ro- vibrational 
and rotational H2 lines (2 x 10 42 erg s " 1 ) is the luminous 
infrared galaxy (LIRG) NGC 6240 (1 Joseph fc Wright I 



[19841: iLutz et al. 112 003: ArmuseLaU|2006|). KeckAdap 
tive Optics and HST NICMOS imaging in near-IR H 2 ro- 
vibrational lines reveal a complex network of hot molec- 
ular gas filaments embedde d in a larger network of Ha 
emission (Ma x et al. Il2005[). This ga laxy has a binary 
active nucleus dKomossa et al. 1 12003) and is thought to 



be a late stage major merger of two gas-rich spiral galax- 
ies. Shock-heating driven by the merger may power the 
st rong H 2 emission spectr um. 

lAppleton et al. I ([20061 ) find high luminosity H2 emis- 
sion from an intergalactic shock between colliding galax- 
ies in Stephan's Quintet. The X-ray emission from the 
shock, though conspicuous, is too weak to power the H2 
lines. Instead, it is likely that the H2 emission is pow- 
ered by kinetic energy dissipat ion through turbulenc e in 
a multi-phase post-shock layer (jGuillard et al. I[2009( ) . In 
this case the H2 emission provides ten times the luminos- 
ity of the next brightest Mid-IR line [Si 11] , demonstrating 
that it is likely the domina nt coolant in this group- wide 
shock (jCluver et al. II2010D . 

H2 ro-vibrational emission from 1000-2000 K molec- 
ular gas is seen from the b rightest central galax- 
ies in some galaxy clusters ([Donahue et al. I [2000). 
Spitzer observations have shown these galaxies also 
display rotational lines from H2 at lower (100-1500 
K) te mperatures ([Egami et al. I 2006; J ohnstone et al. I 
[2001 . The warm molecular gas is often as- 
sociated with extended optical emission line fila- 
ments. Extensive networks of such filaments are 
found in 30-50% of X-ray luminous, rich galaxy clus- 
ters, ass o ciated with the central- d omina n t (cD) galaxy 
(iLvndsl 119701: IHeckman et al. I 119811: iCowie et al. I 



19831 Uohnstone et al. I 119871 : iHeckman et al. I 119891 : 
Grawford et al. II1999D . 



The origin, confinement, and excitation mechanism 
for clu ster emission line filaments are longstanding mys - 
tcrics ( Mi nkowski I [l957t iBurbidge fc Burbidge I |l965). 
The filaments in Perseus A may be formed in-situ by 
the interaction of the radio lobes with the hot intra- 
cluster medium ([Conselice et al. 1 120011 ) . Since the fil- 
aments are dusty and have near-solar metal abundances, 
it is more likely that they are composed of cool gas 
dredged up by the radio jet from the center of the galaxy 
(|Churazov et al. I [20001 iFerland et al~l [20091 ) . Galaxy 
collisions and ram-pressure stripping are two other ways 



to c reate intracluster filaments of atomic and mole cular 
gas (|Kennev et al. 1120081: ISivanandam et al. Il2010f ). 

The gas filaments surrounding Per A are ~ 100 
times more luminous than expected from a cool- 
ing flow and there fore must be powered externally 
(jCowie et al. I 119831 ). The Ha filaments consist of 
extremely narro w threads, most lik ely confined by a 
magnetic field ([Fabian et al. I [2008). Kinetic energy 
dissipation by magnetohydrodynamic (MHD) shocks 
or dissipative MHD waves driven by the radio lobes 
may po wer the atomic and molecular emission line 
spectra (Kent fc Sargent 1 119791: IHeckman et al. 1 119891: 
lHatch et al. I 120051: iGuillard et al. 1 120091) . Molecular 
gas heating by cosm ic ray s has also been suggested 
(|Ferland et all 120081 |2009f), and may help to explain 
some anomalous atomic forbidden emission line ratios. 

1.3. Overview 

We present a large Spitzer spectroscopic survey of 55 
radio galaxies. We find a large number of H2 luminous 
galaxies (17), which belong to and further define the MO- 
HEG class. We use mid-IR and optical forbidden lines, 
aromatic features, and 24 /im continuum emission to in- 
vestigate the gas and dust content and H2 power source 
in these radio-selected MOHEGs. We describe the radio 
galaxy sample, observations, and spectra in Sections 2-4. 
We estimate warm H2 masses using excitation diagrams 
in Section 5. We explore the dust emission properties of 
the sample in Sections 6 and 7. We discuss the H2/PAH 
diagnostic ratio and use this to redefine the MOHEG 
class in Sections 8 and 9. We estimate star formation 
rates and use ionic forbidden line ratios to diagnose AGN 
activity in Sections 10 and 11. We use Chandra obser- 
vations to rule out X-ray heating by the AGN in Section 
12. The host galaxy morphologies and environments are 
discussed in Section 13. Radio jet morphologies are con- 
sidered in Section 14. Finally, we discuss evidence in 
support of radio jet feedback in Section 15.1, cosmic ray 
heating in Section 15.2, and the origin of the molecular 
gas in Section 15.3. 

2. SAMPLE 



We conducted a survey of iFanaroff k Rilev I (Il97l 
class-II radio galaxies an d quasars from the 3CRR cata- 
log w ith the Spitzer IRS ([Ogle, Whvsong. fc Antonucci I 
120061) and a simila r survey of FR I radio galaxies 
( Leipski et al. I [20091 ) . The FR II galaxies were selected 
from the 3CRR catalog to form a redshift and flux- 
limited sample with z < 1 and S„(178 MHz)> 16.4 Jy. 
The FR I galaxies form a representative (yet incomplete) 
sample at z < 0.13, selected from the 3CR catalog with 
a flux limit of 5^(178 MHz)> 15 Jy. 

From the above two samples, we selected all 23 FR 

I galaxies and the subsample of 19 FR II galaxies at 
z < 0.22. We excluded the blazars BL Lac and 3C 
371, which have mid-IR spectra completely dominated by 
synchrotron emission, and the radio-intermediate quasar 
E1821+643. We retrieved the data for 8 additional FR 

II and 5 FR I radio galaxies in the Spitzer archive, for 
a total of 27 FR II and 28 FR I galaxies at z < 0.22. 
The redshift limit of z < 0.22 ensures that the full H2 
S(0)-S(7) line series falls in the Spitzer IRS wave band. 

The radio galaxy sample is presented in Tables 1 and 
2. Alternate names, host galaxy properties, cluster mem- 
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bership, and interaction signatures are listed in Table 3 
for H 2 detected radio galaxies. Throughout this paper 
we assume a cosmology with Hubble constant Hq = 70 
km s _1 Mpc -1 , matter density parameter 51 = 0.3, and 
dark energy density Q\ = 0.7. 

3. OBSERVATIONS 
3.1. Spitzer 

The radio galaxies in Tables 1 and 2 were observed 
with the Spitzer IR Spectrograph (IRS; IHouck et al. I 
12001 . Spectr a were taken with the low-resolution (R ~ 
60-120) modules Short-Low 2 (SL2), Short-Low 1 (SL1), 
Long-Low 2 (LL2), and Long-Low 1 (LL1) over the wave- 
length range 5-36.5 /xm. At the galaxy redshifts of 0.002- 
0.22, l'.'O corresponds to 0.04-3.6 kpc. The slit widths are 
3'.'6, 3'.'7, 10'.'5, and 10'.'7, respectively, corresponding to 
0.15-13 kpc for SL and 0.42-38 kpc for LL. The telescope 
was nodded to place the nucleus of each galaxy at the 
1/3 and 2/3 positions along each slit. 

We began our analysis with the S15-S17 pipeline pro- 
cessed, basic calibrated data sets (BCDs). Off-slit ob- 
servations were subtracted to remove background light. 
Spectra were extracted (using SPICE 2.1) within tapered 
regions matching the Spitzer point-spread function and 
including the first Airy ring (SL2: 7'.'2 at 6 fim, SL1: 14'.'4 
at 12 /im, LL2: 21'.'7 at 16 /im, LL1: 36'.'6 at 27 /xm). We 
tried b oth regular and optim al extraction within these re- 
gions (Na rron et al. Il2007[ ). and selected the extraction 
with higher signal-to-noise for further analysis. For un- 
resolved sources, the continuum and emission line fluxes 
from regular and optimal extractions are in agreement. 

The flux uncertainties produced by the Spitzer pipeline 
are known to be too large for faint sources. The pipeline 
(versions S18.7 and earlier) does not properly adjust the 
uncertainties for 'droop', a systematic, additive effect 
which is proportional to the total photocurrent incident 
upon the detector. We divided flux uncertainties for each 
spectral order by a constant factor of 0.3-2.1, to better re- 
flect the observed noise characteristics. The uncertainty 
scale factor for each order was determined from the ratio 
of the mean pipeline uncertainties to the observed noise 
in emission line-free regions. 

Spectral orders with low flux were adjusted upwards by 
factors of 1.0-2.0 to correct for intcrorder slit loss differ- 
ences. Such differences may arise from pointing inaccura- 
cies, even for unresolved sources. The differences in slit 
widths and wavelength-dependent extraction apertures 
may also lead to mismatches in continuum or line flux 
between orders for spatially resolved sources. Our spec- 
tral order matching will properly correct for this effect in 
the case that the extended continuum and line emission 
are cospatial and have similar spectra, such as might be 
the case for extended star-forming regions. However, the 
emission line equivalent widths and line flux ratios be- 
tween orders could be off by as much as the interorder 
correction factor (< 2.0) in the case of extended nebu- 
lar or H2 emission. This should not be an issue for any 
galaxies except Per A, 3C 293, and Cen A, the only ones 
where the emission lines and PAH features appear to be 
spatially extended along the spectrograph slits. For all 
other sources, the mid-IR (MIR) emission (including the 
H2 emission) is unresolved and resides within the central 
2-40 kpc of the host galaxies. 



We measured the emission line and PAH features of 
all galaxies in a sy stematic way using the spectral fit- 
ting tool PAHFIT (jSmith et al. 1 120071 ). The continuum 
is modeled with a sum of starlight and thermal (black- 
body) emission components, absorbed by a uniform dust 
screen. This is intended only as an empirical description 
of the continuum, not a physical model. The spectrally 
resolved PAH features are fit by Drude profiles. The un- 
resolved H2 and forbidden atomic emission lines are fit 
by Gaussian profiles. Emission line and PAH fluxes from 
the galaxies are presented in Tables 4-8. We also mea- 
sured narrow-band 24 [im fluxes (3.7<7-clipped average 
over 22.5-25.5 /im rest; Tables 1-2) from the IRS spectra 
and used them to estimate mid-IR continuum luminosi- 
ties (Table 9). 

In principle, the continuum for each source should be 
fit with a dust radiative transfer model, but this is out- 
side the scope of the present work. Instead, we tai- 
lored the continuum model in two cases where the sim- 
ple multi-temperature black body plus screen extinction 
model did not work. For Perseus A, we added silicate 
emission features to the continuum. We introduced the 
capability to fit optically thin silicate emission features 
with PAHFIT, computing the silicate emission as black- 
body emission multiplied by an empirical dust extinction 
curve. In the Cygnus A spectrum, 18 /urn silicate ab- 
sorption is unusually weak or absent relative to the 9 /im 
silicate absorption, perhaps because of radiative transfer 
effects. Therefore we r emoved the 18 um silica te absorp- 
tion feature from the iChiar fc Tielensl (j2006l ) Galactic 
center extinction curve in order to model the absorption. 

3.2. Chandra 

In order to test the possibility of H2 heating by X- 
rays from AGN, we compiled unabsorbed 2-10 keV nu- 
clear X-ray luminosities measured by Chandra (Table 9). 
For several sources without published X-ray luminosi- 
ties, we used archival Chandra data to measure them 
ourselves (Table 10). We reprocessed the Chandra data 
using CIAO v4.1.2 with the CALDB v4.1.2 to create a 
new level-2 events file filtered for the grades 0, 2, 3, 4, 
and 6 and with the 0.5" pixel randomization removed. 
To check for periods of high background, we extracted 
light curves for the chip on which the source was placed, 
excluding point sources. Few periods of high background 
were found in any observation, and so we used almost all 
available data in our analysis. 

We extracted the nuclear spectra for each source from a 
small circle (typically of radius 5 pixels or 2.5"), and sam- 
pled background from an adjacent region free from point 
sources. In most cases, the point-like emission from the 
nucleus was easily separated from any extended emission 
components, though there may be significant contami- 
nation from the surrounding hot cluster gas in the case 
of the Zw3146 BCG. We do not consider the possibility 
of H2 heating by X-ray emission from the extended in- 
tracluster medium or X-ray emission from the extended 
radio source. 

We adopted the standard X-ray spectral models for 
radio-loud AGN described by Evans et al. (2006) to fit 
each spectrum. In the case of narrow-line radio galax- 
ies (NLRGs) , the canonical model takes the form of the 
combination of a heavily abs orbed power law and a soft, 
unabsorbed power law that lEvans et all (|2006[ ) inter- 
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prctcd to originate at the base of an unresolved jet. In 
Table 10 we present an observation log and summarize 
the best-fitting models for each source. 

4. H 2 SURVEY RESULTS AND SPECTRA 

We detect two or more H2 emission lines from 17/55 
(31%) of the radio galaxies in our sample (Table 4). We 
require > 2.5a detection of 2 lines at the systemic red- 
shift. A typical 2a H2 S(l) detection limit for the sample 
is 8 x ltT 16 erg s" 1 cm" 2 (8 x 1CT 19 W m" 2 ), but this 
may vary with redshift and H2 line contrast relative to 
the continuum. 

H2 line and PAH fluxes for galaxies with only one or no 
H 2 lines detected are presented in Table 6. Upper limits 
to their H2 luminosities are estimated as Lh2 < 1.6L[H2 
0-0 S(l)+ S(3)], where the correction factor 1.6 accounts 
for the S(0) + S(2) line flux and is the mean value of [H2 
S(0) + S(l) + S(2) + S(3)]/[H 2 S(l) + S(3)] for galaxies 
where all 4 lines were detected. We use this rather than 
the sum of 4 upper limits for S(l) through S(4), which 
would be an overestimate. 

The Spitzer IRS spectra of the H2 detections are pre- 
sented in Figures 1-7, along with archival HST optical im- 
ages for most sour ces. Some of these spectra were prev i- 
ously presen t ed bylQgle. Whvsong. &: Antonucci I (|2006f l , 
lOgle et al. I (|2007| ). and iLeipski et al. I (|2009f ). but we 
show them here as well for complete ness. The Cen A 
spectr um was originally presented by iWeedman et al.~l 
(2005). H2 pure-rotational lines, forbidden lines of O, 
Ne, Si, S, Ar, CI, and Fe, and polycyclic aromatic hydro- 
carbon (PAH) features are seen in many of the spectra 
(Tables 4-8). Galaxies with large 7.7 /im PAH equiva- 
lent width are presented in Figure 1 and those with small 
7.7 fxm PAH equivalent width in Figures 2-3. The mid- 
IR luminous radio galaxies Per A, Cyg A, 3C 433, and 
3C 436, which contain dust-obscured quasars, are pre- 
sented in Figure 4. The H2 emission is difficult to see on 
top of the continuum emission from the galactic nuclei 
in some spectra. However, subtraction of a continuum 
model clearly reveals the H2, PAH, and other emission 
lines (Per A, Cyg A, and Cen A Figs. 5-7). 

5. H 2 LUMINOSITIES, TEMPERATURES, MASSES, AND 
SURFACE DENSITIES 

In order to compute the mass of warm molecular gas, it 
is first necessary to know the H2 excitation mechanism. 
The low critical densities (H2-H2 n cr it < 10 3 cm _3 at 
T > 100 K) of the H 2 0-0 S(0) and S(l) emissi on lines 
(|Le Bourlot et al. Ill999t iWrathmall et al. 1120071 ) ensure 
that their upper levels will be collisionally excited in al- 
most any environment where warm H2 is found. Non- 
thermal excitatio n by UV photons or cosmic-ray sec- 
ondary electrons (|Ferland et al. I [20081 ) may contribute 
to the higher rotational and ro- vibrational levels, but 
are inefficient at exciting the lower pure-rotational levels. 
Therefore the lowest H 2 energy levels, which dominate 
the warm H 2 mass budget, must be in thermal equilib- 
rium. 

Assuming thermal equilibrium in gas at a single tem- 
perature T, the H2 luminosity is determined by the H2 
mass Affj2 and T according to 

^(Ha) = f^H 2 (T) 
where the mean luminosity of a single H2 molecule, 
summed over all transitions u — > I is 



'(Ha) = E»,i hv ul A ul f u (T) 
The occupation number f u (T) of each level is determined 
by the Boltzmann factor and partition function. For ex- 
ample, a molecular gas mass of Mg = MH 2 /lO 9 Af and 
a temperature of 200 K gives 

L(H 2 ) = 1.1 x 10 42 erg s" 1 M 9 
For multiple temperature components, the luminosity of 
each component can be computed separately. 

For all radio galaxies with two or more H2 emission 
lines detected, we construct excitation diagrams to de- 
termine the temperature distribution and mass of H2. 
For the limiting case of an H2 source size of 3'.' 7 x 3'.' 7 
uniformly filling the width of the SL slit, we plot the col- 
umn density of the upper level divided by the statistical 
weight for each transition (Figs. 8-10). 

We simultaneously fit the data with 3 temperature 
components (Table 11), the minimum number necessary 
to fit the emission line fluxes and parameterize what is 
likely a continuous temperature distribution. We quote 
la single parameter uncertainties for the H2 tempera- 
tures and masses. Upper limits are included in the fits 
as 2a detections with la errorbars. This allows us to 
give conservative upper limits for the lowest tempera- 
ture component for sources where the S(0) line is not 
detected, without affecting much the other components. 
Note that the lowest temperature component, which is 
constrained by the S(0) and S(l) lines dominates the to- 
tal H2 mass. We give an upper limit for the warmest 
temperature component mass in 3C 338, where the S(6) 
and S(7) lines are not detected. 

We assume that the H 2 level populations and ratios 
of ortho/para H2 are in local thermodynamic equilib- 
rium (LTE) and adjust the statistical weights of the or- 
tho transitions iteratively to obtain a self-consistent fit. 
There is no evidence for non-equilibrium ortho/para ra- 
tios, which would manifest as suppressed or enhanced 
emission from odd-J relative to even- J transitions. The 
ortho/para ratios may depend upon the detailed ther- 
mal history of the gas, and can take much longer than 
the gas cooli ng time (~ 10 4 yr) to reach their equilib- 
rium values (iWilgcnbu s et al. 1 12000). However, for the 
warm (T > 100 K) H 2 observed by Spitzer, the equilib- 
rium ortho/para ratio takes on a narrow range of values 
(1.6-3.0) such that the equilibrium assumption has only 
minor consequences for our model fits. For example, as- 
suming an ortho/para ratio of 3.0 rather than 2.1 de- 
creases the temperature of the coolest H2 component in 
3C 326 from 110 to 100 K, and raises the mass by only 
27% to 2.8 ±0.6 x 10 9 M ? . 

To determine the sensitivity of our model fits to the 
number of temperature components, we tried both 2- 
componcnt and 3-component fits to all of the lines in the 
3C 326 N excitation diagram. The 2-component fit gives 
a mass of 8.7 ± 3.2 x 10 8 M© at 140 K (x 2 /DOF = 36.6), 
compared to 2.2± 1.2 x 10 9 M Q at 110 K (x 2 /DOF = 1.8) 
for the 3-component fit. The 3-componcnt fit has much 
smaller residuals and is clearly better. The additonal 
temperature component increases the estimated H2 mass 
by a factor of 2.5 in this case. Which lines are detected 
and used in the fit can also change the estimated H2 
mass. If we throw out the S(5), S(6), and S(7) lines and 
refit the 2 component model, we find 1.6 ± 0.7 x 10 9 M Q 
at 120 K (x 2 /DOF = 1.9). 
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We measure radio galaxy warm H2 masses from 8 x 
10 6 Af Q up to an amazing 2 x 10 10 M Q , the bulk at tem- 
peratures of 100-200 K (Tables 11 and 12). Such large H 2 
masses imply the wholesale heating of an entire galaxy's 
worth of molecular gas in some cases. Smaller masses 
(3 x 10 4 -8 x 10 7 M Q ) of warmer H 2 at 200-1500 K are also 
found in all galaxies where H2 is detected. Comparable 
luminosity is emitted by each temperature component in 
a given galaxy, though most of the mass resides in the 
coolest component (Table 11). In all cases where there 
are several high-level lines detected, the H2 temperature 
increases with upper-level excitation temperature. This 
produces the characteristic curves in Figures 8-10, which 
have less negative slopes at higher energies. This behav- 
ior is also s een in the excitation dia gram for the Perseus A 
filaments (jJohnstone et al."ll2007t ). This is expected for 
any astrophysical situation wher e there is a range of te m- 
peratures and densities present (jFerland et al. 1 12008). 

The total warm H2 column densities and column den- 
sity lower limits range from < 4 x 10 20 — 1 x 10 23 cm~ 2 
(3 - 8OOM pc~ 2 ). Only lower limits can be inferred for 
spatially unresolved H2 emission. This spans the range 
observed for normal star-for ming galaxies at the low en d 
to ULIRGs at the high end (jGracia-Carpio et al. II2007D . 

CO emission has been detected by past observations 
(lEvans et al. lll999HLim et al |[200l'(Salome fc Combes I 
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Saripalli fc MackF g007; 
JNesvadba et al. I 12010 ; 
Ocana Flacquer et al. I I2010|) of several radio MO- 
HEGs (Table 12). The mass of cold (10 - 50 K) H 2 
estimated assuming the standard Galactic CO to H 2 
conversion factor of 4.6 M /(K km s _1 pc 2 ) ranges 
from 4 x 10 6 — 2 x 10 10 M Q . However, we caution 
that the applicability of the Galactic conversion factor 
to the extreme conditions in radio MOHEGs has not 
been established. This caveat notwithstanding, the 
inferred ratio of warm (T > 100 K) to cold H2 in these 
galaxies ranges from < 0.18-2.1. The galaxies 3C 218, 
272.1, 326, 424, and 433 are particularly remarkable in 
having warm/cold H2 ratios of > 0.5. A large fraction 
of the H2 in these galaxies is disturbed and heated to 
temperatures > 100 K. 

6. 24 fiM DUST EMISSION 

The formation and shielding of H2 requires the pres- 
ence of dust, which may reveal it self through IR contin - 
uum and PAH emission features (|Guillard et al. lf20Toh . 
However, there must be a sufficient UV radiation field 
to heat the very small dust grains (VSGs) and excite 
the PAH stretching or bending modes in order to ob- 
serve mid-IR dust signatures. A UV photon field may be 
sustained by 1) star formation, 2) AGN emission, or 3) 
thermal bremsstrahlung and line em ission from radiative 
shocks (|Dopita fc Sutherland II 1996h ■ It must be kept in 
mind that the H2, VSG, and PAH emission regions are 
spatially unresolved and may or may not be coincident. 
Also, different types of sources can contribute in different 
ratios to each of these spectral features. 

First we consider the ratio of H2 t o 24 /im dust emissio n 
(Fig. 11) for star-forming galaxies ([Roussel et al. Il2007l ). 
There is a trend (with large scatter) for Li^^jL ^A to 
decrease with L24. ULIRGs (Higdo n et al. Il2006l ) have 
smaller L(H2)/£24 than nearby spirals and dwarfs, by a 



factor of ~ 5. This is similar to the known trend of de- 
creasing L(P AH)/L 2 4 with I/24 observed for starb ursts 
and ULIRGs (|Desai et al. I [20061: iTran et al. |[200ll and 
Fig. 12). It is debated whether this latter trend is in- 
trinsic to star-forming regions or whether it reflects an 
increasing AGN contribution to the continuum at higher 
luminosity. In either case, the tendency of H2 and PAH 
emission to decrease together relative to L24 may sug- 
gest that they are both powered by star fo rmation in 
starbu rsts and ULIRGs. On the other hand, Zakamska 
([201CA finds that H 2 emission is less obscured than PAH 
emission in ULIRGs, indicating a different spatial distri- 
bution. 

The radio galaxies in our sample have up to 20 times 
larger ratio of L(H2)/L24 than do normal star-forming 
galaxies (Fig. 11). In some cases (e.g., 3C 326 and 
Stephan's Quintet), > 10% of the bolometric IR lu - 
minosity comes out in the H2 lines (O gle et al. I [20071. 
This dramatically indicates that the H2 emission is not 
powered by star formation, and that H2 emission is an 
important and perhaps primary coolant for the ISM in 
these sources. Strong H2 emission is often an indicator 
of magnetic (C-ty pe) shocks in Galacti c sources such as 
Herbig-Haro jets (iNeufeld et al. Il2006f) an d supernovae 
(jCesarskv et al. lll999HNeufeld et al. II2007I ). A deficit of 
24 /im continuum relative to H2 emission indicates the 
importance of non-radiative (e.g., shock) heating. 

The SINGS AGNs also have enhanced L(U 2 )/L 2i 
([Roussel et al. 1120071 ) . though they are not as well sepa- 
rated from the star-forming galaxy locus as radio galax- 
ies (Fig. 11). Non-radiative heating of H2 appears to 
be important in these nearby, radio-quiet AGNs. A sim- 
ilar enhancement in L(H 2 )/Z/ 2 4 is seen for some of the 
iKaneda et all ([20081 ) 'dusty elli pticals', which are RSA 
ellipt icals detected by IRAS ([Goudfrooii fc de Jong I 
119951 ). Inter estingly, 12/14 of th e H 2 detected ellip- 
ticals in the IKaneda et al. I ([20081 ) sample have radio 
sources detected in the NVSS or SUMSS survey-7 com- 
pact sources and 5 FR Is. Radio jet activity could be 
responsible for exciting the excess H2 emission in these 
sources. 

7. PAH EMISSION 

PAH emission features are detected in all of the ra- 
dio MOHEGs except Cyg A (Table 5) and in half of 
the H2 nondetected radio galaxies (Table 6). How- 
ever, the ratio of 7.7 /im PAH to 24 /im contin- 
uum emission L(PAH7.7)/L24 is generally lower by 
a factor of 10-100 compared to normal star-forming 
galaxies from the SINGS sample (Fig. 12). The 
mid-IR continuum of most radio galaxies and ra- 
dio MOHEGs is therefore not star-formation domi- 
nated. In many of these galaxies, the continuum has 
been shown to be AGN dominated, including contri- 
butions from AGN -heated dust and synchrotron emis- 
sion from the iet (lOgle. Whvsong. fc Antonuccil 120061: 
Cleaiy et aL~ll2007t iTadhunter et al. Il2007t lHaas et al. I 
20081: iLeipski et al. (l2009f T 

The large scatter that we see in i(PAH7.7) /L24 (2 dex) 
primarily reflects varying star formation vs. AGN frac- 
tional contributions to the spectra. Additional scatter is 
introduced by variation in the 7.7/11.3 /im PAH band 
ratio (Fig. 13). Low L(PAH7 '.7) / L24 is also observed 
in the IKaneda et al. 1 (|2008|) dusty ellipticals, indicating 
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that the warm dust continuum in these sources does not 
come from star formation. 

The median PAH(6.2 /xm)/PAH(7.7 /zm) ratio is 0.25 
for radio MOHEGs, which is indistinguishable from 
the median ratio (0.28) for normal star-forming SINGS 
galaxies (Fig. 13). The median PAH(7.7 /^m)/PAH(11.3 
fim) ratio of 1.8 is significantly lower in radio MO- 
HEGs than the median value of 4.1 measured in star- 
forming SINGS galaxies and with in the range of va l- 
ues reported for SINGS AGNs (jSmith et al. I l2007h . 
The iKaneda et al. I (|2008l ) dusty elliptical galaxies have 
a similar median 7 .7 to 11.3 /im band ratio of 1.4. 
iSmith et al. I (2007) present an extensive discussion of 
possible interpretations of the observed differences in 7.7 
to 11.3 fim band ratio between AGNs and star-forming 
galaxies and among AGNs. 

The impact of the PAH ionization state and size dis- 
tribution on P AH emission spectra has been quanti- 
fied in models (iBakes et al. I [200ll: iDraine fc Lil[200ll: 
iFlaeev et al. 112006ft . In these models, the 7.7 to 11.3 ^m 
band ratio is mainly sensitive to the ionization state and 
is lower for neutral PAHs than for cations. This inter- 
pretation is supported by observations of photodissocia- 
tion regions (PDRs) which show that the 7.7 to 11.3 /jm 
band ratio decreases wi th increasing depth into the cloud 
(|Rapacioli et al. I [2005') . The ionization states of PAHs 
depend on the PAH ionization parameter Guv x /n e) 
where G V y is the 912 < A < 2400 A UV field en- 
ergy density r elative to the Solar Neighborhood value 
(jHabing l[l96 8). and T e and n e the temperature and elec- 
tron density. The electron density is that of the neutral, 
molecular, and atomic gas where PAHs are present. 

The interpretation of Galactic data leads to equal frac- 
tions of neutral and cation PAHs in the diffuse Galactic 
ISM for a mean ionization parameter ~ 1000 K 5 cm 3 
and a 7.7 to 11.3 /xm band ratio (m easured with PAH- 
FIT) of 2.9±0.2 (|Flagev et al. $ 2006). The median value 
of the 7.7 to 11.3 fxm band ratio of radio MOHEGs (1.8) 
is smaller than the value for the Galactic diffuse inter- 
stellar medium. This implies a smaller PAH ionization 
parameter, which could be accounted for by lower Guv 
in galaxies with a low star formation rate, or by larger 
n e from increased ionization by cosmic rays or X-rays. 

In the Stephan's Quintet intergalactic shock, the ob- 
served PAH emission is consistent with the mass of warm 
H2 gas heated by a radiation field with Guv = 1-4, as- 
suming a Galactic dust-to -gas mass ratio and a galac- 
tic dust size distribution (jGuillard et al. I [2010D . This 
shows that PAHs can survive in dense regions of the 
interstellar medium which are shocked at velocities < 
100 km s _1 . PAHs could potentially be destroyed by 
faster shocks, which ma y be present in the diffuse ISM 
(jMicelotta et al. Il2010al) . In the hot 10 6 - 10 7 K X-ray 
emitting gas phase, PAHs would mostly be destroyed by 
electron collisi ons and would not sur vive more than a 
few 10 3 years (jMicelotta et al. 112010b] ). PAHs are likely 
protected from electron sputtering in the denser, colder 
phases. 

We similarly expect to see PAH emission from the large 
H2 reservoir in radio MOHEGs, excited by the ambient 
UV radiation field. PAH destruction by X-rays may pos- 
sibly reduce PAH abundances in the vicinity of the AGN 
(jVoit II 19911 ). However, the four radio MOHEGs with sig- 



nificant star formation (discussed in §10-3C 31, 293, 315, 
and Hydra A) and correspondingly strong UV radiation 
field show nearly normal 7.7 and 11.3 /urn PAH equiv- 
alent widths and PAH band ratios, demonstrating that 
PAH abundances have not been depleted by the AGNs 
on a galactic scale. The remaining radio MOHEGs show 
relatively weak PAH emission, consistent with a weaker 
UV field. 

8. H 2 TO PAH EMISSION RATIO 

We plot the ratio of H 2 luminosity in the 0-0 S(0)-S(3) 
lines over L(PAH7.7) versus L 2 4 (Fig. 14 ). Normal star- 
formi ng galaxies from the SINGS survey (|Roussel et al. I 
12001 have a median ratio L(H 2 )/L(PAH7.7) = 0.014, 
which appears to be independent of L24 over four decades 
(L 24 = 10 6 - 10 10 Lq). The mean L(H 2 )/L(PAH7.7) 
has a similar but slightly greater val ue for ULIRGs 
(jHigdon et al. Il2006t iDesai et al. 112006ft . 

In contrast, all of the H 2 detected radio galaxies have 
very large L(H 2 )/L(PAH7.7) = 0.03 - 4 or greater. The 
mo st extreme source s, with L(H 2 )/Lpah > 3 are 3C 326 
N (|Ogle et al. II2007D and 3C 424. Other ga laxies of dif- 
feren t types, including the LIRG NGC 6240 dLutz et al. I 
12003ft . the Zw 3146 BCG (|Egami et al. 112006ft. and the 
Steph an's Quintet intergalactic shock (jAppleton et al. I 
2006) are also unusual in having large L(H 2 )/Lpah ra_ 
tios. Strikingly, all of t he SINGS Seyferts an d all but 4 
of the SINGs LINERs (jRoussel et al. 1 120071 ) also stand 
out as having relatively large L(H 2 )/Lpah = 0.02 — 0.2, 
though not so large as the most extreme radio galax- 
ies. Such large ratios indicate shock, cosmic ray, or 
X-ray heating rather than UV heating: of the molecu- 
lar g as phase (jRigopoulou et al.~l 120021 : iRoussel et al. I 
120071 ). We explore the various possible heating mech- 
anisms for H 2 in radio galaxies below. 

Alternatively, it might be argued that a high ratio 
L(H 2 )/L(PAH7.7) is caused by a low abundance of PAHs 
in the warm molecular gas rather than low UV field. A 
UV field of Guv ~ 4 would be required explain the 11.3 
/im PAH flux from 3C 326 N if it is associated with the 
observed warm H 2 and CO (Nesvad ba et al. 1 12010T ) . as- 
suming a normal H 2 /CO conversion factor (Table 12). 
We use archival GALEX UV data to independently con- 
strain the intensity of the radiation field and PAH abun- 
dance. The NUV AB magnitude is 21.47, which trans- 
lates into a surface brightness of 1.59 x 10 -6 W m~ 2 sr _1 , 
assuming emission from star-forming regions uniformly 
covering a 5 arcsec (8 kpc) diameter area. This corre- 
sponds to a mean UV radiation field Guv = 12. If any 
of the observed UV radiation comes from the AGN, the 
Guv value contributing to PAH excitation would be less. 
A comparison of the two Guv estimates gives an upper 
limit on the PAH depletion of ~ 3. This is far from suffi- 
cient to explain the extreme ratio L(H 2 )/L(PAH7.7) > 4 
in 3C 326 N, which is 300 times the mean value in normal 
star-forming galaxies. 

9. MOHEG DEFINITION 

Since we do not have far-IR measurements for most of 
the sources in ou r sample, we do no t use our original MO- 
HEG definition (jQgle et al. II2007D . We could re-express 
that definition as L(H 2 )/L 2 4 > 2 x 10~ 2 (dashed line in 
Fig. 11), excluding all normal star-forming galaxies from 
the MOHEG class. The observed H 2 to continuum con- 
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trast ratio for our H2 detected radio galaxy sample ranges 
from L(H 2 )/L 2 4 — 6 x 10~ 4 — 0.4, which means that only 
9/17 of these galaxies would be considered MOHEGs by 
such a definition. Inspection of Figure 11 suggests that 
defining MOHEGs in terms of a minimum i(H 2 )/i24 
contrast ratio is problematic. For one, this ratio de- 
creases with Z/24 for normal star-forming galaxies, so a 
line with negative slope in that figure would do a bet- 
ter job at separating MOHEGs. However, it would not 
completely account for the variation in AGN continuum 
contribution at 24 /im, which introduces scatter into this 
relation. 

The much cleaner separation between MOHEGs and 
normal star-forming galaxies in Figure 14 compared to 
Figure 11 motivates us to redefine MOHEGs as galaxies 
with L(H 2 )/L(PAH7.7) > 0.04. This particular cutoff 
was chosen to exclude all non-AGN SINGs galaxies from 
the MOHEG class. With our new definiti on, 16/17 of the 
H 2 de tected radio galaxies, 10/18 of the iKaneda et al.~l 
(2008) dusty ellipticals, NGC 6240, and the Zw 3146 
BCG are all MOHEGs. We call radio galaxies that 
are MOHEGs 'radio MOHEGs' to distinguish them from 
MOHEG samples selected in other wavebands. The radio 
galaxy 3C 31, with J L(H 2 )/L(PAH7.7) = 0.03 falls just 
below the MOHEG cutoff as a result of the relatively 
large contribution of star formation and strong 7.7 fim 
PAH emission its mid-IR spectrum (Fig. 1). However, 
for simplicity and since it has a larger L(H 2 )/L(PAH7.7) 
ratio than any SINGs galaxy of similar luminosity Z 2 4, 
we will also call this source a radio MOHEG. It is not 
possible with the current data to reliably classify the H 2 
nondetected radio galaxies with respect to the MOHEG 
criterion. Therefore, the 31% MOHEG fraction in our 
radio galaxy sample should be considered a lower limit. 

The success of the L(H 2 )/L(PAH7.7) diagnostic ratio 
in distinguishing between AGN and star-formation dom- 
inated spectra reflects the different molecular gas heat- 
ing mechanisms in the two classes. Furthermore, the 
L(H 2 )/L(PAH7.7) value at a given L24 might be used 
to determine the ratio of AGN to star-forming contribu- 
tions to these two emission features. If the H 2 and 7.7 
fim PAH emission have different spatial extent, this ra- 
tio and classification as a MOHEG might depend upon 
the aperture used to extract the spectrum. There is a 
similar issue with AGN optical classification, since AGN 
fraction tends to increase with decreasing aperture size. 

10. STAR FORMATION 

In general, 3C radio galaxies tend to have less lumi- 
nous PAH emission and lower inferred star-formation 
rates than other AGN sa mples, including PG quasars 
and 2MASS red quasars ([Shi et al. I 120071 ) . Whether 
this is separ ately true for those with and withou t hid- 
den quasars (|Ogle. Whysong. fe Antonucci 1120061 ) needs 
to be investigated. For only 4/17 radio MOHEGs do the 
spectra (Fig. 1) indicate a large star formation fraction 
(3C 31, 218, 293, and 315), with L(PAH7.7)/£ 2 4 > 0.1 
(Fig. 12). Based on template fitting, the star formation 
contributi on in 3C 31, 218 and 293 ranges from 87-100% 
at 15 fim (|Leipski et al. Il2009r i. Using the same proce- 
dure, we find that 3C 315 has a 40-50% star formation 
contribution at 15 fim. 

We use the observed 7.7 /im PAH luminosity to es- 
timate star formation rates (SFRs, Table 13) for ra- 



dio M OHEGs using the prescription of iRoussel et al. I 
(2001). We convert the 11.3 fim PAH luminosity to 7.7 
fim PAH luminosi ty, dividing by the mean 11.3/7.7 /im 
PAH ratio of 0.26 ([Smith et al. 1120071 ) to make an equiv- 
alent SFR estimate for sources where the 7.7 fim PAH is 
not detected. Note that the SFRs are integrated over 
the Spitzer IRS slit, which does not cover the whole 
galaxy for most of the observations (Figs. 1-4). The 
SFRs for MOHEGs with detected 7.7 /jm PAH emission 
range from 4 x 10~ 3 Af Q yr" 1 for 3C 272.1 at the low 
end to 2.6M yr -1 for 3C 436 on the high end. Since 
the 11.3/7.7 fim PAH ratio is systematically enhanced 
compared to normal star- forming galaxies (Fig. 13), the 
11.3 /im PAH star-formation rates should be considered 
upper limits and used with caution. 

11. AGN ACTIVITY 

Most (13/17) of the radio MOHEGs a r e mid-IR weak 
B ([Ogle. Whvsong. fe Antonucci I 120061 : iLeipski et aT~l 
120091) . with vL v {2^fim) < 2 x 1O 1O L (Table 9). Most 
are also optically classified as low-ionization galaxies 
(LIGs, Table 1) based on their low [Om] 5007/[On] 
3727 rati os. Using the opt i cal em ission line measure- 
ments of iButtiglione et al. I (120091) .we co nstruct BPT 
diagnostic diagrams ( Baldwin et al. |[l981l ). which show 
that most radi o MOH EGs are LINERs according to the 
iKewlev et al. I (|2006| ) line ratio criteria (Fig. 15). We 
can not say whether this tendency of radio MOHEGs 
to host MIR-weak, LINER AGNs is intrinsic or rather 
represents a bias against detecting low-equivalent width 
H 2 emission lines and PAH features against the stronger 
continuum in the MIR-luminous HIGs and BLRGs. 

The mid-IR forbidden line ratios of radio MOHEG 
LINERs (Fig. 16) are similar to those of other LIN- 
ERs observed by Spitzer, with [Ne m]/[Ne il]< 2.0, [O 
ivl/ [Ne nl < 1-0, JS ml / JNe 11R 1.0 and [S ml /[Si n < 



0.5 (iSturm et al. 1120061: iDale et al. ll2lMlDudik et al. I 
I2009TI . The low [S III] /[Si 11] ratio in particul ar is charac- 
terist ic of AGN X-ray dissociation regions ([Dale et al. I 
2006), in contrast to the higher ratio typical of star- 
forming galaxies. The [Ne m]/[Ne 11] and [S m]/[Si 11] ra- 
tios of radio MOHEGs are however considerably greater 
than what are seen i n the Stephan's Quintet main shock 
(jCluver et al. 112010ft . This may indicate that radio MO- 
HEG forbidden line emission is associated primarily with 
the AGN rather than the shocked H 2 emission region. 
On the other hand, we can not rule out the possibility of 
a significant contribution from fast J-type shocks (with 
similar line ratios to LINERs) to the optical and MIR 
forbid den emission lines, (e.g., 3C 326. iNesvadba et al. I 
l201Ch . 

The three radio MOHEGs with Seyfert-like, high- 
ionization optical spectra (Cyg A, 3C 433, and 436) 
all have [Ne m]/[Ne n]> 2.0 and [O iv]/[Ne 11] > 
2.0. High ionization [Ne v] and [Ne vi] lines are de- 
tected in Cyg A and 3C 433, but only the [Ne vi] 
7.65 itm line is detected in 3C 436 (Table 8). All 
three host powerful, hidden quasars with vL^(24/im) > 
2 x 1O 1O L , and deep silicate troughs (T9.7 > 0.5) in- 



1 T he wavelength used in the IQg lc. Whvsong, & Ant onucci I 
pOM ) definition of MIR-weak {vL v (15Atm) < 8 X 10 43 erg s-1 
= 2 X 10 10 -Lq) is different than the one we employ here, because 
of the need to accomodate a different redshift range. 
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dicate obscuration of the AGN by a dditional cold dust 
(jQgle. Whvsong. fc Antonucci I [2006). Polarized, scat- 
tered broa d emission lines ar e seen in the optical spectra 
of Cyg A (jQgle et al. II1997T) and 3C 433 (Ogle et al., in 
preparation), confirming the presence of hidden quasar 
nuclei. The [Ne v | and [Ne Vij lines are detected at 2-3cr 
inCen A (see also lWeedman et al. 1120051 : iMirabel et al. I 
1999), which may contain a hidden Seyfert nu cleus 
(jWhvsong fc Antonucci! I2004t iLeipski et al. 1120091) . 

12. INSUFFICIENCY OF AGN X-RAY HEATING AS AN H 2 
POWER SOURCE 

The lack of a high-luminosity, high-ionization AGN in 
most radio MOHEGs indicates weak accretion onto the 
central supermassive black hole. Such low luminosity 
AGNs are insufficient to power the H2 emission via X- 
rays. 

The process of H2 heating by X-rays in X-ray disso- 
ciation regions (XDRs) is described by iMalonev et al. I 
(fl996l) . In their XDR models, 30-40% of the absorbed 
X-ray flux goes into gas heating via photoelectrons. The 
atomic photoelectric cross section energy dependence is 
E~ 8 / 3 , so that photoionizations in the energy range 1- 
30 keV are important. For a standard AGN X-ray 
power law spectrum with frequency dependence j/~ ' 7 , 
L(l - 30keV)/L(2 - lOkeV) = 2.3. 

For a characteristic observed H2 temperature of 200 
K (corresponding to an ionization parameter log£ e ff = 
—2.6) the cooling by H2 rotati onal lines in XDR mo dels is 
~ 2% of the total gas cooling (jMalonev et al. Ifl996l Figs. 
3a, 5). At this temperature, the ratio of the first four 
rotational lines to the total rotational line luminosity is 
L(H 2 0-0 S(0)-S(3))/i(H 2 )= 0.58. Combining the above 
factors, we estimate a maximum H2 to X-ray luminosity 
ratio of L(H 2 0-0 S(0)-S(3))/i x (2-10 keV)< 0.01. This 
ratio is conservative, since it assumes that all of the X-ray 
flux from the AGN is absorbed by the XDR. 

We compare unabsorbed nuclear X-ray (2-10 keV) 
luminosities measured by Chandra (Table 9) with the 
summed H2 S(0)-S(3) line luminosities measured by 
Spitzer (Fig. 17). Two sources lack X-ray nuclear flux 
measurements (3C 310 and 424). We find L(H 2 )/i x (2- 
10 keV)> 0.01 in all measured radio MOHEGs and 
SINGS AGNs, except Cyg A. The warm H 2 in Cyg A 
would have to intercept and absorb > 30% of the X- 
ray luminosity from the AGN to produce the observed 
mid-IR H2 line luminosity. For all other sources, the ob- 
served rotational L(H2)/Lx(2-10 keV) exceeds the maxi- 
mum XDR H2 to X-ray luminosity ratio estimated above. 
The large scatter (nearly 5 dex) in the observed ratio 
L(H 2 )/L X (2-10 keV)= 0.002 - 100 also argues strongly 
against X-ray heating as the primary driver of H2 emis- 
sion in radio galaxies or SINGS AGNs. 

13. HOST GALAXIES AND ENVIRONS 

Many radio MOHEGs (Ta ble 3 and Figs. 1-4) a ppear 
to have a galactic st ellar disk (iDonzelli et al. 112007ft , dust 
disk, or dust lanes (|De Koff et al. 1120001 ) . Stellar or dust 
disks may indicate a past gas-rich merger, which would 
help to explain the presence of large quantities of H 2 . 
Most of the other galaxies are peculiar or distorted in 
their stellar isophotes. Only 3C 386 appears to have a 
regular elliptical host morphology. 

There is evidence for patchy or organized dust absorp- 



tion in 12/17 radio MOHEGs (Table 3). This is im- 
portant because H2 forms most easily in the presence of 
dust. Dust masses m easured from HST optical images 
(|De Koff et al. I I2000D are typically < 1% of the warm 
H2 mass measured with Spitzer (Table 12). However, 
this method likely underestimates the total dust mass 
associated with the warm H2 and cold CO. The hosts 
of 3C 31 and 270 have striking dust disks and Cen A 
has thick dust lanes from an edge-on warped disk. Dust 
lanes, patchy dust, or dust filaments are also seen in Per 
A, 3C 272.1, 293, 315, 317, 338, Cyg A, 3C 433, and 
436. The only galaxies without obvious dust absorption 
are Hydra A, 3C 310, 326 N, 386, and 424. However, 
dust absorption could be missed in the HST snapshots 
because of a smooth dust distribution or low galaxy in- 
clination. 

Interestingly, 14/17 radio MOHEGs belong to close 
galaxy pairs, groups or clusters (Table 3). Of these, 3 
reside in pairs, 3 in groups, 1 in a poor X-ray detected 
cluster, and 7 in X-ray luminous clusters. Four of the 
clusters (hosting Per A, Hydra A, 3C 317, and 3C 338) 
are so-called cool-core clusters with ~ 3 keV gas temper- 
atures and high central X-ray surface brightness. The 
cooling flow phenomenon may be a way to deliver large 
quantities of gas to these radio MOHEGs. FR I radio 
galaxi es tend to reside in rich environments at low red - 
shift (jLongair fc Seldner I Usmk iHeckman et al. 1 119851) . 
while low-z FR II radio galaxies prefer the field. Be- 
cause radio MOHEGs are found in both environments, 
cluster membership is apparently not a prerequisite for 
the H2 emission phenomenon. Gravitational interactions 
with nearby companions may play a significant role in 
delivering or driving the molecular gas into the center of 
radio MOHEGs in the field or in clusters. A statistical 
comparison of the frequency of physical galaxy pairs in 
the H2 detected and nondetected radio galaxy subsam- 
ples is outside the scope of this paper. 

There are tidal distortions indicating strong interac- 
tions in at least 4/7 galaxy pairs, groups, and poor clus- 
ters: 3C 293, 310, 326, and 433. The 3C 31 and 3C 315 
pairs do not show strong distortion, but the halos of close 
companions appear to overlap the primary galaxy. The 
Virg o cluster galaxy 3 C 270 (NGC 4261) has a faint tidal 
tail (jTal et al. II2009D . The 4 cool-core cluster cDs and 3C 
433 have large, diffuse halos which appear to envelop two 
or more galaxies in the cluster centers. The three galax- 
ies that do not reside in pairs, groups or clusters are 3C 
386, 436, and Cen A. The warped, dusty molecular disk 
in Cen A may be the result of an elliptical/spiral merger 
(|Baade fc Minkowski I [l95l IQuillen efahl 1200(1 . The 
host of 3C 436 shows signs of a recent merger, includ- 
ing a second n ucleus and other irreg ularities in its sur- 
face brightness (jMcLure et al. I fl9991 . The Virgo cluster 
galaxy 3C 272.1 (M84), the group galaxy 3C 424, and the 
isolated galaxy 3C 386 show no obvious signs of recent 
interaction. 

14. RADIO SOURCE MORPHOLOGIES 

The radio MOHEGs in our sample cover a large range 
of radio luminosity i/L l/ (178MHz) = 10 39 - 10 44 erg s" 1 
(Table 9 and Fig. 18). They display a motley collec- 
tion of radio morphologies (Tables 1 and 2). There are 
6 FR lis and 11 FR Is. The fraction of FR II MOHEGs 
is comparable to the fraction of FR I MOHEGs in the 
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sample (22 ± 9% vs. 39 ± 12 %). However, there could 
potentially be a bias against detecting H2 in FR lis with 
hidden quasars because the emission lines might be over- 
whelmed by copious warm dust emission. 

Further classification of radio morphologies is some- 
what subjective but may offer clues to the origin of radio 
MOHEGs. Looking at their morphology sub-classifica- 
tions (Table 1), 5/11 FR Is are classical twin jet (TJ) 
sources, 4/11 are fat doubles (FD), 1/11 is a double- 
double (DD), and 1/11 has a compact symmetric core 
(CSC). Of the FR lis, only 3/6 are classical doubles, 1/6 
is a fat double (FD), and 2/6 are X-shaped (X). There 
is a relatively large fraction of distorted radio sources 
compared to the rest of the sample, which only has 2/38 
FD and 0/38 X-shaped sources. The two X-shaped ra- 
dio sources (3C 433 and 3C 315) and the very peculiar 
source 3C 338 may be distorted by the motions of their 
interacting hosts or have jets deflected by their compan- 
ions. 

The radio MOHEGs also include a high fraction 
(3/17) of so-called restarting double-double (DD) or CSC 
sources, compared to 1/38 of the H2 nondetected radio 
galaxies. The radio map of 3C 424 shows a brighter in- 
ner set of lobes surrounded by a fainter set of lobes that 
appears to origina te from a previous spurt of activity 
(Blac k~e"t al. Ill992h . High resolution maps of 3C 293 and 
315 reveal bright, inner, compact symmetric double radio 
lobes roughly aligned with the larger scale radio lobes. 
These may represent new outbursts of radio activity or 
interaction of the small scale radio source with the host 
galaxy ISM. The radio cores also appear relatively bright 
in 3C 310, and 386, perhaps also indicating rejuvenated 
jets, or else high Doppler beaming factor. The restart- 
ing jet phenomenon may be connected to a new supply 
of fuel provided by the large mass of H2 in these galax- 
ies. Interaction between the restarting jet and the host 
galaxy ISM may play an important role in heating the 
H 2 . 

15. DISCUSSION 
15.1. Radio-Jet Feedback 

The high incidence of MOHEGs in our radio-selected 
sample may indicate that the radio jets or lobes power 
the H2 emission. However, there is at best a weak cor- 
relation between radio MOHEG H2 luminosity and 178 
MHz radio lobe power Lns (Fig. 18a). The linear Pear- 
son correlation coefficient is R = 0.79, which reduces 
to R = 0.35 (< 17% probability of a chance correla- 
tion) when removing the partial correlation of luminos- 
ity with distance. For MOHEGs, there is a large range 
L(H2) I Lns — 0.004 — 6 in the ratio of H 2 luminosity 
to low-frequency radio lobe synchrotron luminosity. The 
upper limits for H 2 nondetected radio galaxies have sim- 
ilar scatter, with a lower median ratio than the radio 
MOHEGs in this plot. The available mass of H 2 must 
be a limiting factor, and is likely an important source of 
scatter in the H 2 luminosity. 

Scatter in the 178 MHz radio luminosity may be caused 
by environmental density variations or spectral aging of 
the radio lobes. Size and time scale differences between 
the radio and H 2 emission regions may also be impor- 
tant. For most sources, the spatially unresolved (kpc 
scale) H 2 emission region is much smaller than the ra- 



dio source, which extends far beyond the host galaxy (50 
kpc-1 Mpc). The extended radio emission may not give a 
reliable indicator of the jet power dissipated at kpc scale. 
However, while the H2 cooling time scale is only ~ 10 4 
yr, the dissipation timescale of the molecular gas kinetic 
energy may be even greater than the typical radio source 
lifetime of - 10 7 vr (jNesvadba et al. Il2010t) . Therefore, 
variations in jet power over the lifetime of the radio jet 
may not cause significant variations in the H2 luminosity. 

There is no significant correlation between L(H2) and 
the 5 GHz radio core power L5 measured at arsecond 
scales (Fig. 18b). This is not surprising since relativistic 
beaming can boost or de-boost the core flux by a large 
factor. Differences in radio spectral shape may also be 
large at this high frequency. Neither is there a correlation 
ofL(H 2 )/£(PAH7.7)with J Li 78 norL 5 (Fig. 18c,d). This 
ratio and the fraction of warm/cold H2 may depend on 
a number of factors, including the spatial distribution of 
H2 relative to the radio source, the pressure of the hot 
ISM, and the detailed jet-ISM interaction history. 

It is instructive to compare radio MOHEGs to dusty 
elliptical MOHEGs selected by IRAS (jKaneda et al. I 
120081 ). We estimate 178 MHz luminosities for 12 ellip- 
ticals by interpolating or extrapolating 100 — 10 3 MHz 
radio fluxes measured by the NVSS, SUMSS, or FIRST 
surveys, or fluxes available in NED (Fig. 18a, c). These 
galaxies are up to 100 times fainter than radio MO- 
HEGs at 178 MHz, and extend the MOHEG L(H 2 ) vs. 
Li78 parameter space. The faintest radio sources with 
Lira < 10 38 erg s^ 1 (in NGC 2974, NGC 3894, IC 3370, 
and NGC 5044) lie well above the L(H 2 ) vs L 17S locus 
established by radio MOHEGs, with L(H 2 )/Li 78 > 50. 
We suggest that these weak radio sources may be radio 
galaxies that are currently in between strong radio jet 
outbursts. This hypothesis can be tested by searching 
for low- frequency radio relics. The H2 emission may be 
powered by residual kinetic energy injected into the ISM 
by previous radio jet activity. The H2 luminosity floor of 
L(H 2 )~ 10 39 erg s _1 may be set by the ratio of the H2 
kinetic energy dissipation timescale and the mean time 
between radio outbursts. Alternatively, this low level of 
H2 emission could be powered by gas accretion or X-ray 
emission from a weak AGN. It will be important to test 
the latter possibility with additional Chandra observa- 
tions. 

Several of the radio MOHEGs that reside in X-ray 
bright clusters show evidence that the radio lobes have 
displaced the intracluster medium (ICM) in X-ray cavi- 
ties. For these sources, it is possible to estimate the 4pV 
energy required to excavate the cavities, and divide by 
the buoyancy time scale to estimate the averag e jet power 
deposited into the ICM (jRaffertv et al. 1120061 ). These jet 
cavity powers fall in the range 10 42 — 10 45 erg s _1 and 
the ratio of H2 luminosity to jet cavity power falls in the 
range 1 x 10~ 4 - 3 x 10~ 3 (Table 14). This demonstrates 
that the jet power is more than sufficient to produce the 
observed H2 luminosity, modulo an unknown conversion 
efficiency. 

One potential difficulty is efficiently coupling the jet 
power to the host galaxy ISM after the radio lobes have 
broken out of the host galaxy and are doing most of t heir 
work on the intergalactic medium. lOgle et al. I (|2007|) ar- 
gued that only 7 x 10~ 4 of the radio jet power in 3C 326 
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could be dissipated in the host galaxy ISM, based on the 
small ratio of radio core to Mpc-scale radio lobe flux. 
This in turn could only power < 10% of the observed 
H2 luminosity, for an estimated jet kinetic luminosity of 
~ 10 44 erg s _1 . Given the uncertainties in estimating jet 
kinetic luminosity and time variability, the inner radio 
jet could easily be ten times more powerful than this es- 
timate. In addition, the radio core power is probably not 
a reliable tracer of jet kinetic power dissipation and may 
be strongly affected by Doppler boosting or de-boosting, 
as indicated by the large scatter in Figure 18b. Turn- 
ing this problem around, the large H2 luminosity from 
the center of 3C 326 N may indicate that a significant 
fraction of the radio jet power is dissipated inside the 
host galaxy, even though the radio lobes have long since 
escaped. 

On the other hand, the compact sy mmetric radio 
lobes i n the radio M O HEG 3C 293 (lEvans et al I 
1991 lEmonts et al. I l2005t IQgle et al. I 120081 
Papadopoulos et al. I |2010T) appear to be directly 
interacting with the host ISM. It will be important to 
measure the pressure in the hot interstellar medium 
with Chandra to determine if dissipation of jet kinetic 
energy can indeed power the H2 emission in this and 
other kpc-size sources (e.g., 3C 315). 

High-velocity neutral and ionized outflows have been 
observed in the radio MOHEG 3C 293, and there is 
morphological and kinematic evi dence that these out- 
flows are driven by t he radio jet (jMorganti et al. I [20031 : 
lEmonts et al. I 120051 ) . Similar H 1 and ionized out- 
flows have been found in several other low-redshift com- 
pact and extende d radio sources (|Morganti et al. I [20051 : 
iHolt et al. 112008ft . with outflow speeds of up to 2000 km 
s _1 and mass-loss rates of up to 50 Mq yr" 1 indicat- 
ing a strong inter action between t h e rad io jet and the 
host galaxy ISM. INesvadba et al. I ([2010D recently dis- 
covered a high velocity (~ 1000 km s _1 ) neutral (Na 
I) outflow in 3C 326, the most extreme radio MOHEG. 
Even more spectacular jet driven ionized outflows, which 
may be responsible for ejecting a large fraction of the host 
galaxy ISM, are observ ed in high-redshift radio galaxies 
(|Nesvadba etaT~l[2008[ ). 

The connection between jet driven outflows and H2 
emission in rad i o MO HEGs is studied in detail by 
INesvadba et""aT~l ((2010ft . Since the molecular hydrogen 
must have high density and a low volume filling factor, it 
is likely that the cocoon of hot, shocked ISM surro unding 
the jet (jScheuer 111971 iBegelman fc Cioffi l[l989h serves 
as an intermediary between the jet and the H2. The ra- 
dio jet heats the diffuse ISM, which expands in a hot fast 
wind and drives shocks into ambient or entrained molec- 
ular clouds, powering the observed warm H2 emission. 
For H 2 at densities of 10 3 -10 4 cm~ 3 , the observed range 
of temper atures is characterist ic of 5-20 km s" 1 shock 
velocities (|Guillard et al. II2009D . 

The relative velocity and shear between the molecular 
clouds and hot wind may ultimately po wer the H2 emis- 
sion through turbulent mi xing layers (jNesvadba et alH 
l2010HGuillard et al. Il2009ft . The hydrodynamical stabil- 
ity of the H 2 emitting clouds needs to be investigated, as 
they may be ablated by the hot wind and may eventually 
dissolve completely. However, there is evidence that the 
atomic and H2 filaments in Per A are held together by 



strong magnetic fields (iFabian et al. I [20081) . Moreover, 
the fact that H2 is the primary coolant only in magnetic 
(C-type) shocks strongly suggests that the magnetic field 
plays an important role in MOHEGs. 

15.2. Cosmic Ray Heating and Pressure 

Another potentially important source of heating for 
warm H2 in radio MOHEGs may be ionizing particles 
in the form of r elatiy istic particles from the radio lobes. 
iFerland et al. I (|2008ft have proposed cosmic rays or ther- 
mal particles from the surrounding hot gas as the dom- 
inant heat source powering the warm H2 emission fila- 
ments in the Perseus A (3C 84) cool core cluster. An 
enhanced cosmic ray density has also been proposed to 
account for the high tempera ture of molecular clouds i n 
the Milky Way nuclear disk (|Yusef-Zadeh et al. I [20071) . 
We discuss this possibility for H 2 luminous radio galax- 
ies. 

The critica l density of the S(0) and S(l) lines is suf- 
ficien tly low (|Le Bourlot et al. 1119991 iWrathmall et al. I 
120071 < 10 3 cm" 3 at T > 100 K) that we can assume 
that the corresponding excited levels are thermalized by 
collisions. Thus, the line emission is cooling the gas and 
must be balanced by gas heating. In our radio MOHEG 
sample, the mean of total line emission per warm H2 
molecule (in the lowest temperature component, Table 
11) is observed to be 7 x 10~ 24 ergs -1 . 

The molecular gas heating by cosmic rays is 4 eV per 
ioniz ation in the MEUDON gas code (|Le Petit et al. I 
2006). In CLOUDY, the heating efficiency increases with 
the gas ionization fraction; the heating per ionization is 
7 eV for a ionization fraction of 10~ 4 characteristic of 
the diffuse interste llar medium where carbon is ionized 
(|Shaw et al. I [2005). The H 2 line cooling is balanced by 
cosmic ray heating for an ionization rate per hydrogen 
( ~ 7 x 10~ 13 s _1 , a value 3.5 x 10 4 times h igher than 
the standard Galactic rate (|Shaw et al. 1 12005). For such 
a value, cosmic rays are the main destruction path of H2 
molecules and the molecular gas fraction depends on the 
ionization rate to g as density ratio C/uh- Model calcu- 
lations presented bv lFerland et al.1 ([2008) show that for 
£ ~ 7 x 10~ 13 s _1 the gas is molecular for njj > 5 x 10 4 
cm~ 3 (see their figure 2). Note that for such high den- 
sities H2 rotational states are thermalized up to J > 5, 
and the higher J lines up to at least S(3) are also cooling 
lines. 

Combining the lower limit on the gas density and the 
warm H2 temperature (~ 200 K) inferred from the S(l) 
to S(0) line ratio we get a warm molecular gas pressure 
~ 10 7 Kcm~ 3 which is a factor 1400 higher than the 
mean cos mic ray pressure of 7.2 x 1 3 K cm -3 in the Milky 
Way disk (|Boulares & Coxlll990f) . For a cosmic ray en- 
ergy density 3.5 x 10 4 times the Galactic mean value, 
the cosmic ray pressure would be 25 times larger than 
the thermal gas pressure (assuming that the cosmic ray 
pressure scales linearly with the ionization rate). Such 
a large difference in pressure would seem unsustainable, 
however it could in principle be supported by magnetic 
pressure in t he compressed magnet ic field within molec- 
ular clouds (Pelle grini et al. 1 120091 ) . Therefore, heating 
by cosmic rays from the radio lobes or hot gas in the ra- 
dio jet cocoon may present an alternative or additional 
pathway to heat H 2 by radio-jet feedback. 
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15.3. Gas Supply from Cooling Flows and Galaxy 
Collisions 

An important unanswered question is the origin of 
large quantities of molecular gas in radio MOHEGs. 
Copious amounts (10 9 — 1O 11 M ) of cold molecular 
gas discovered in ~ 45% of BCGs in X-ray selected 
cool-core clusters (jEdge 1120011 ISalome fc Combes ll200l 
may be deposited by cooling flows. H2 emission from 
warm molec ular gas also appears to be common in cool- 
core BCGs dEgami et al. l l2006t Uohnstone et al. 1 120071: 
iDe Messieres et al. 112009ft . Since a large fraction of these 
BCGs host radio sources, the H 2 emission could poten- 
tially be excited by radio jet feedback. The gas cool- 
ing rates in the radio MOHEGs t hat reside in cool core 
clusters are 10 - 200 My yr" 1 ([Peterson et al. 1 120031 : 
iHicks fc Mushotzky 1 12005ft , enough to supply the ob- 
served molecular gas masses in < 1 Gyr. 

Gas-rich galaxy collisions or tidal interactions with 
nearby galaxies may drive gas into the center of radi o 
galaxies in less rich environments ( Ogl e et al. 1 12007ft . 
This is an attractive possibility, considering that most 
radio MOHEGs are mergers or have nearby compan- 
ions, and several display prominent tidal distortions. In 
comparison, most field elliptical galaxies and 3C radio 
galaxies in the nearby universe are poor in molecular 
gas (< 10 9 M^. iCombes et al. I [20071: iLim et al. I l200l 
lOcana Flacquer et al. 1 12010ft . Galaxy collision driven 
inflows and jet induced outflows may constitute an- 
other, intermittent type of AGN feedback loop (e.g., 
lAppleton et~aLll2002ft . 

While cooling flows or galaxy collisions may deliver the 
molecular gas in radio MOHEGs, it seems unlikely that 
gas accretion or galaxy collisions contribute significantly 
to heating the H2. Unlike the galaxy- collision powered 
H 2 e mission from Stephan's Quintet (jAppleton et al. I 
2006), the bulk of the H 2 emission is concentrated inside 
the central ~ 10 kpc of radio galaxies. We are pursuing 
higher resolution spectroscopy of mid-IR and near IR H 2 
emission lines in order to determine the dynamical state 
of the warm molecular gas. 

16. CONCLUSIONS 

(1.) We detect high luminosity H 2 emission lines from 
17/55 (31% of) radio galaxies at z < 0.22 observed with 
Spitzer IRS at low resolution. 

(2.) We redefine molecular hydrogen emission galax- 
ies (MOHEGs) to be galaxies with L(H 2 )/L(PAH7.7) > 
0.04. This puts 16/17 of the H 2 -detected radio galaxies 
(radio MOHEGs) in this rapidly growing, new class of 
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galaxies. 

(3.) Large H 2 /PAH ratios indicate that the H 2 emis- 
sion from MOHEGs is most likely powered by shocks. We 
can not completely rule out a contribution from cosmic 
ray heating, though it would require cosmic ray pressures 
that exceed the gas pressure by a large factor. 

(4.) The star formation rates in radio MOHEGs range 
from very low to moderate (4 x 10 -3 — 3M Q yr _1 ), and 
the 24 /jm continuum is dominated by the AGN in 13/17 
of these galaxies. 

(5.) Most MOHEG AGNs in our sample are mid- 
IR weak and have low-ionization forbidden emission line 
spectra. 

(6.) AGN X-ray emission measured by Chandra is not 
luminous enough to power the H 2 emission in MOHEGs. 

(7.) All but one of the radio MOHEGs in our sam- 
ple reside in merger remnants, interacting pairs, groups, 
or clusters, supporting the hypothesis that galaxy inter- 
actions and mergers may deliver molecular gas or drive 
existing molecular gas into the central few kpc of radio 
galaxy hosts. 

(8.) Radio jet driven outflows of hot gas may drive 
shocks into the molecular gas, powering the H 2 emission 
in radio MOHEGs. The jet power measured for radio 
MOHEGs in cool core clusters is more than sufficient to 
power the H 2 emission. Higher spatial resolution of the 
H 2 emission region and its kinematics will give further 
insights into the radio jet feedback mechanism. 
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Fig. 1. — Spitzer IRS low-resolution spectra of 3C 31, 218 (= Hydra A), 293, and 315. These sources have strong PAH emission and 9.7 
fim silicate absorption. Spectra are arranged in order of increasing equivalent width of the H2 pure-rotational lines. Images to the right 
are from HST WFPC2 (3C 31, 293 and 315) and the UK Schmidt telescope (3C 218, from NED). The host of 3C 31 resides in the Arp 
331 group. Hydra A resides in the cool core cluster A 780. The box at the lower left corner of each image shows the size of the 3'.' 7 X ICK'0 
spectral extraction region at 10 fim. 
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Fig. 2.— Spitzer IRS low-resolution spectra of 3C 270 (= NGC 4261), 3C 272.1 (= M 84), 3C 310, and 3C 317. These sources have weak 
7.7/im PAH emission and weak-moderate 9.7 fim silicate absorption. Note the extremely large PAH 11.3/PAH 7.7 and PAH 17/PAH 7.7 
ratios and steeply rising continuum at > 24 fim in 3C 272.1. Radio galaxies 3C 270 and 272.1 reside in the Virgo Cluster, 3C 310 in a 
poor cluster, and 3C 317 in the A 2052 cool-core cluster. The box at the lower left corner of each image shows the size of the 3"7 X lO'.'O 
spectral extraction region at 10 (tm. 
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Fig. 3. — Spitzer IRS low-resolution spectra of 3C 326, 338, 386, and 424. These sources have we ak PAH emission and no silicate 
absorption. HST WFPC2 and ACS images are to the right. The 3C 326 pair appears to be interacting (Ogle ct al. 2007). Radio galaxy 
3C 338 resides in the A 2199 cool core cluster, while 3C 386 is an isolated elliptical galaxy. The host of 3C 424 appears (centered on the 
smaller PC chip) in a rich group of galaxies. The box at the lower left corner of each image shows the size of the 3'/7 X 10''0 spectral 
extraction region at 10 fj,m. 



Molecular Hydrogen in Radio Galaxies 



17 



15 



I 

5 



H 2 7 6 £(5J 5(4) 5(3) 5(2) S(l) 
[Fe ll] ' [AMI] Wlll]'[S IV] [Ne ll] [Ne III]' [S ill] 



S(0) [SMI] 
[0 IV] ' [5 Iff] 



PAHS.Z 7.7 B.6 11,3 12.6 
I II II 



3C 436 








■ 

■ & 


1 


■ . 


i - 






<=>„ 






Fig. 4. — Spitzer IRS low-resolution spectra of 3C 436, Cyg A, 3C 84 (= Per A), and 3C 433. Silicate emission is seen at 10 (im and 18 
fim in the Per A spectrum, and silicate absorption in 3C 436, 433, and Cyg A. HST WFPC2 images are at the right. Per A is interacting 
with an in-falling dusty companion seen in the foreground. The disky, dusty host of 3C 433 resides in an interacting group of 3 galaxies. 
The box at the lower left corner of each image shows the size of the 3'.'7 X lO'.'O spectral extraction region at 10 fim. See Figs. 5 and 6 for 
spectral fits of 3C 84 and Cyg A. 
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Fig. 5. — Spitzer IRS low-resolution spectrum of 3C 84 (= Per A), with best fit model and model components. Silicate emission at 10 (im 
and 18 /im is fit by a blackbody times an empirical AGN silicate dust emissivity curve. The bottom panel shows the continuum-subtracted 
PAH and emission line spectrum, including H2 pure-rotational lines. 
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Fig. 6. — Spitzer IRS low-resolution spectrum of 3C 405 (= Cyg A), with best fit model and model components. A custom extinction 
curve (dotted line), lacking the 18 fim silicate feature, was required to fit the continuum. The bottom panel shows the continuum-subtracted 
emission line spectrum, including H2 pure-rotational lines. 
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Fig. 7. — Spitzer IRS low-resolution spectrum of Cen A (= NGC 5128), with best fit model and model components. The bottom panel 
shows the continuum-subtracted PAH and emission line spectrum, including H2 pure-rotational lines. 
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Fig. 8. — H2 excitation diagrams. Upper level column density divided by statistical weight (for H2 uniformly filling the SL slit width) or 
2<t upper limit is plotted against upper level energy E/k^(K) for each H2 0-0 S(J) pure-rotational emission line. Model H2 temperature 
components (dotted lines) are listed in Table 11. 
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Fig. 9. — H2 excitation diagrams. Upper level column density divided by statistical weight (for H2 uniformly filling the SL slit width) or 
2<r upper limit is plotted against upper level energy E/k^(K) for each H2 0-0 S(J) pure-rotational emission line. Model H2 temperature 
components (dotted lines) are listed in Table 11. 



Molecular Hydrogen in Radio Galaxies 



23 




Fig. 10. — H2 excitation diagram for Cen A. Upper level column density divided by statistical weight (for H2 uniformly filling the SL 
slit width) or 2cr upper limit is plotted against upper level energy E/ki(K) for each H2 0-0 S(J) pure-rotational emission line. Model H2 
temperature components (dotted lines) are listed in Table 11. 
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Fig. 11. — Ratio of H2 luminosity summed over the 0-0 S(0)-S(3) pure rotational lines to L24 = vL u (2A fim , rest) lum inosity. Radio 
galax ies (this work) are compared to SINGS st ar-forming galaxies, Seyferts, and LINERs (Rousscl ct al. 2007), ULIRGs (Higdon ct al. 
12009) . and the subset of IKaneda et al. I H2008I) dusty ellipticals that have 24 /an Spitzer MIPS fluxes IITemi et al. Il2007f). Radio galaxy 
H2 upper limits (from T ables 2 and 6) are pl otted as downward pointing arrows. The LIRG NGC 6240 ( Armus et al. 2006), the Zw 3146 
brightest cluster galaxy (jEgami ct al. 2006), and the Stephan's Quintet intergalactic shock (SQ shock sub-region. ICluver et al~l 120101 ) 
are plotted for comparison. The large values of HYI2) / L24 > 0.02 in many radio galaxies indicate the importance of non-radiative heating 
(e.g., shocks). 
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Fig. 12. — PAH to vL v (2A fim, rest) luminosity ratio, for the 7.7 and 11.3 /an PAH features. (See Fig. 11 for symbol legend; H2 
nondetected radio galaxies are plotted as open triangles.) Soli d line corresponds to 1:1 ratio. Rad io galaxies (this work) are compared to 
SINGS star-forming galaxies, Seyferts, and LINERs flRoussel et al~ll2007l: ISmith et al~1l2007l ') ■ the lKaneda et al. I {2008) dusty ellipticals, 
and the Stephan's Quintet intergalactic shock (SQ shock sub-region, ICluver et al. 1 120101, Most radio galaxies have L(PAH7.7)/L24 < 0.1, 
indicating a very small star formation contribution and dominant AGN contribution to the mid-IR continuum. 
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Fig. 13. — Left: PAH emission feature flux ratios (see Fig. 11 f or symbol legend). H2-d etected radio g a laxies (triangles, this work) are 
compared to SINGS star-forming galaxies, Seyferts, and LINERs (Smith et al. 2007), the Kancda ct al. (2008) dusty ellipticals, and the 
Stephan's Quintet intergalactic shock (SQ shock sub-region, Cluver et al. 2010;). Four radio galaxies with neither 6.2 fim nor 7.7 fim PAH 
detected are not plotted. The median PAH(7.7 £tm)/PAH(11.3 fim) ratio is lower in radio galaxies than in star-forming galaxies, while 
the median PAH(6.2 /jm)/PAH(7.7 /im) ratio is similar. This indicates lower PAH ionization caused by a weaker interstellar UV radiation 
field or increased g as-phase ionization by X-rays or cosmic rays. Right: PAH 7.7 fim/PAH 11.3 fim vs. [Ne ni]/[Ne n] line flux ratio after 
ISmith ct al. (200^). One radio galaxy with no PAH features detected (3C 405) is not plotted. The large [Ne lll]/[Ne II] ratio relative to 
star-forming galaxies may be attributed to AGN activity. 
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Fig. 14. — Ratio of H2 luminosity summed over the 0-0 S(0)-S(3) lines to 7.7 (im PAH luminosity vs. vL v (2i /^m, rest) continuum 
luminosity. This ratio indicates the relative importance of mechanical heating and star formation power. All but one of the H2 -detected 
radio galaxies from our sample stand out above normal star-forming galaxies from the SINGS survey {Rousscl et al. 2007) and are 
MOHEGs, with L(H 2 )/L(PAH7.7) > 0.04. Radio galaxy H 2 upper limits with detected 7.7 fira PA H emission (from Tab les 2 and 6) arc 
plo tted as downward poi nting arrows. Many of the LINERs and Seyferts from the SINGS sample (Rousscl et al. 200 ?]) and several o f 
the lKaneda et al. I J20 08) dusty ellipticals are MOHEGs. Other MOHEGs from the literature, includ ing NGC 6240 (lArmus et~al~ll2009) , 
Zw 3146 I IEeami et al. II2006T ). and the Stephan's Quintet intergalactic shock (SQ shock sub-region, ICluver et aT~l l2010), are plotted for 
comparison. 
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Fig. 15. — BPT (Baldwin ct al. 1981) AGN diagnostic diagrams, a) The radio MOHEGs (filled triangles) have optical emission line 
ratios consistent with pure AGN, rather than star-forming (H II) or com posite galaxies. Optical line fluxes for these galaxies and other 
2 < 0.3 3C radio galaxies (small dots) are from Buttiglionc ct al. (2009). b) and c) Nearly all radio MOHEGs on these plots (except 3C 
433 and 436) are classified as LINERs according to the dividing lines of lKewfev et al. I i 20061 ). d) The density-sensitive [S II] doublet ratio 
falls within the normal range for radio galaxies. 
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comparison in the first panel, c) Dale et al. (2006) MIR forbidden line diagnostic diagram (see Fig. 14 for symbol legend). The line ratios 
of radi o MOHEGs with [S i 11] measured are consistent with LINERs and Seyferts, but segregate from star- forming galaxies in t he SINGS 
survey (Smith et al. 2007). The Stephan's Quintet intergalactic shock (red pentagon, SQ shock sub-region, Cluvcr ct al. 2010]) has lower 
ionization than radio MOHEGs. 
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Fig. 18. — a,b) H2 luminosity (summed over 0-0 S(0)-S(3)) vs. radio luminosity at 178 MHz and arcsecond-scale radio core luminosity at 
5 GHz (http://3crr.cxtragalactic.info/cgi/database). Radio MOHEGs (filled triangles), H2 nondetected radio galaxies (downward arrows), 
and Kancda ct al. (2U08) dusty ellipticals (open circles) are shown. Dashed lines correspond to 1:1 ratio. c,d) H2 over 7.7 fim PAH 
luminosity ratio vs. radio luminosity at 178 MHz and arcsecond-scale radio core luminosity at 5 GHz. 
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TABLE 1 

H 2 Detected Radio Galaxies (Radio MOHEGs) 



Source 


Morph. a 


Sp.T. b 




z 


F 2 A C 


S178 d 


Ref. 


3C 31 


LTJ 


LIG 





.01701 


19.9±0.3 


18.3 


1 


3C 84 


I,FD 


LIG 





.01756 


3020. ±30. 


66.8 


1 


3C 218 


I,TJ 


LIG 





.05488 


9.1±0.2 


225.7 


4 


3C 270 


I,TJ 


LIG 





.007465 


42.4±0.4 


53.4 


2 


3C 272.1 


I,TJ 


LIG 





.003536 


29.4±0.2 


21.1 


1 


3C 293 


I,CSC 


LIG 





.04503 


34.3±1.0 


13.8 


1 


3C 317 


I,FD 


LIG 





.03446 


4.5±0.1 


53.4 


3 


3C 338 


Ip,FD 


LIG 





.03035 


4.4±0.3 


51.1 


3 


3C 386 


I,FD 


LIG 





.01689 


2.7±0.4 


26.1 


1 


3C 424 


I,DD 


LIG 





.1270 


2.2±0.1 


15.9 


3 


Ccn A 


I,TJ 


LIG 





.001825 


4950.±110. 


1005. 


1 


3C 310 


II.FD 


LIG 





.0538 


1.9±0.2 


60.1 


1 


3C 315 


II,X,CSC 


HIG 





.1083 


2.5±0.1 


19.4 


1 


3C 326 N 


II 


LIG 





.0895 


0.53±0.12 


22.2 


1 


3C 405 


II 


HIG 





.05608 


800. ±20. 


9483. 


3 


3C 433 


II,X 


HIG 





.1016 


154.± 3. 


61.3 


1 


3C 436 


II 


HIG 





.2145 


4.0±0.1 


19.4 


1 



a Radio morphologies: I or II = Fanaroff & Riley ( 1974) type (p = peculiar), 

TJ = twin jet, FD = fat double, CSC = compact symmetric core, DD = 

double-double, X = 'X'-shaped, SJ = single (one-sided) jet. 

b Optical spectral type: LIG = low ionization narrow line galaxy, HIG = high 

ionization narrow line galaxy, BLRG = broad line radio galaxy. 

c Rest 24 /im flux (mjy) measured with Spitzer IRS. A 3.7<r-clipped average 

over 22.5-25.5 fim was used, removing any large noise spikes or contributions 

from the [Ne v] 24.31 /an line. 

d The 178 MHz flux density (Jy) on lBaars et~al~l (fl977f> sca le. 

e Refe renc es for 178 MHz flux density . (l) ILaing et al. I H1983T) f2')lKuhr et al I 
flTMlT ) (3) IKellermann et al. I pMifi (4) Estimated from IKuhr et alH (19811 ) 
160 MHz flux density. 



Ogle et al. 



TABLE 2 

H2 NONDETECTED RADIO GALAXIES 



oourcG 


l\/T(~lT*T-lVl a 

ivlOl JJIl. 


Cr, T b 
op. 1 . 


z 


-F24 


7K d 
01/0 


Rcf c 


3C 15 


1 


LIG 


0.0730 


12.4±0.3 


16.9 


2 


3C 29 


I,TJ 


LIG 


0.045031 


11.2±0.3 


16.2 


2 


3C 66B 


I,TJ 


LIG 


0.02126 


6.6±0.1 


26.8 


1 


3C 76.1 


I,FD 


LIG 


0.03249 


2.6±0.2 


13.3 


1 


3C 83.1 


I 


LIG 


0.02514 


4.9±0.1 


26.6 


1 


3C 120 


I,SJ 


BLRG 


0.03301 


573. ±5. 


7.4 


2 


3C 129 


I,TJ 


LIG 


0.0208 


7.5±0.2 


51.1 


3 


3C 189 


I,TJ 


LIG 


0.04284 


6.8±0.1 


7.4 


2 


3C 264 


I,TJ 


LIG 


0.02172 


15.2±0.3 


28.3 


1 


3C 274 


I,TJ 


LIG 


0.004360 


75.±1. 


1145. 


1 


3C 318.1 


I 


LIG 


0.04531 


0.9±0.2 


12.3 


3 


3C 348 


I,TJ 


LIG 


0.1540 


0.6±0.2 


387. 


2 


3C 403.1 


I 


LIG 


0.0554 


< 0.2 


14.7 


3 


3C 449 


I,TJ 


LIG 


0.017085 


0.75±0.11 


12.5 


1 


3C 465 


I,TJ 


LIG 


0.03022 


4.3±0.2 


41.2 


1 


NGC 6251 


IP 


LIG 


0.02471 


42. ±2. 


10.9 


1 


IC 4296 


Ip,DD 


LIG 


0.01247 


13.7±0.2 


16.8 


2 


3C 17 


lip 


BLRG 


0.219685 


12.6±0.4 


22.2 


2 


3C 28 


11 


LIG 


0.1953 


< 0.3 


17.8 


1 


3C 33 


11 


HIG 


0.0597 


97.1±0.8 


59.3 


1 


3C 61.1 


11 


HIG 


0.1878 


7.2±0.3 


34.0 


1 


3C 111 


11 


BLRG 


0.0485 


168. ±1. 


70.4 


3 


3C 123 


lip 


LIG 


0.2177 


8.2±0.2 


206.0 


1 


3C 192 


11 


HIG 


0.0597 


9.8±0.4 


23.0 


1 


3C 219 


11 


BLRG 


0.1744 


ll.liO.l 


44.9 


1 


3C 234 


11 


HIG 


0.1848 


289. ±1. 


34.2 


1 


3C 319 


11 


LIG 


0.1920 


0.35±0.10 


16.7 


1 


3C 321 


II 


HIG 


0.0961 


360.±10. 


14.7 


1 


3C 346 


II,FD 


HIG 


0.1620 


9.0±0.1 


11.9 


1 


3C 381 


II 


HIG 


0.1605 


44.2±0.5 


18.1 


1 


3C 382 


11 


BLRG 


0.0579 


91. ±2. 


21.7 


1 


3C 388 


11 


LIG 


0.0917 


1.7±0.1 


26.8 


1 


3C 390.3 


11 


BLRG 


0.0561 


231. ±1. 


51.8 


1 


3C 401 


11 


LIG 


0.2011 


1.5±0.1 


22.8 


1 


3C 445 


11 


BLRG 


0.0562 


243. ±2. 


25.2 


2 


3C 452 


11 


HIG 


0.0811 


57.6±0.3 


59.3 


1 


3C 459 


11 


HIG 


0.2199 


103. ±3. 


30.8 


2 


Pic A 


11 


BLRG 


0.03506 


129.±2. 


411. 


4 



a Radio morphology (see Table 1 for key). 
b Optical spectral type (see Table 1 for key). 

c Rest 24 fim flux (mjy) or 2<r upper limit measured with Spitzer IRS. A 
3.7tr-clipped average over 22.5-25.5 (im was used, removing any large noise 
spikes or contributions from the [Ne v] 24.31 /im line. 

d The 178 MHz flux density (Jy) on lBaars etaT~l JT977I) scal e. 

e References for 1 78 M Hz flux density. ( 1 ) ILaing et al~l JT983T) . (2) 
Kuhr et al."| lUMil l (3) IKellermann et al. I 11969) (4) Estimated from 
Kuhr et al.l 119811) 160 MHz flux density. 
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TABLE 3 

Radio MOHEG Host, Environment, and Interactions 



Source 


Name 


NGC 


Host a 


Dust Morph. b 


Environment 


Interaction 11 


3C 31 




383 


Ep 


disk 


Arp 331 group 


close companion 


3C 84 


Per A 


1275 


cD+D 


complex 


Perseus CC cluster 


infalling spiral 


3C 218 


Hyd A 




cD 


A 780 CC cluster 


close companion 


3C 270 




4261 


Ep 


disk 


Virgo cluster 


tail 


3C 272.1 


M 84 


4374 


Ep 


lanes 


Virgo cluster 




oO 2yo 






on 

bO 


lanes 


pair 


bridge, tail 


3C 310 






Ep+D 




poor cluster 


bridge 


3C 315 






SO 


lane 


pair 


close companion 


3C 317 






cD+D 


filaments 


A2052 CC cluster 


close companion 


3C 326 N 






Ep+D 




pair 


bridge, tail 


3C 338 




6166 


cD 


filaments 


A2199 CC cluster 


nest 


3C 386 






E 




single 




3C 405 


Cyg A 




Ep 


complex 


Cyg A cluster 


merger 


3C 424 






Ep 




group 




3C 433 






SO 


complex 


group 


bridge, nest 


3C 436 






Ep 


lane 


single 


merger 


Cen A 




5128 


Ep 


lanes 


single 


merger 



a Host morphology. E= elliptical, Ep= peculiar elliptical, cD = cluster dominant, SO = lenticu- 
lar. Th e tag '+D' indicates a significant exponential disk component in the near-IR radial light 
profile UDonzelli et al. Il2007l 1. 
b Dust morphology, after JDe Koff et al. Il2000ft . 

c Cluster, group, or pair membership. CC cluster = cool X-ray core cluster. 

d Qualitative assessment of tidal interactions including close companion, tail, bridge, or merger. 

'Nest' indicates multiple galaxies surrounded by a common envelope. 
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TABLE 4 

Radio MOHEG H 2 0-0 S(J) Emission Lines 



Source H 2 S(0) H 2 S(l) H 2 S(2) H 2 S(3) H 2 S(4) H 2 S(5) H 2 S(6) H 2 S(7) 
28.22 fim 17.03 /im 12.28 9.66 jim 8.03 fim 6.91 fim 6.11 (jm 5.55 A»m 



3C 31 
3C 84 
3C 218 
3C 270 
3C 272.1 
3C 293 
3C 310 
3C 315 
3C 317 
3C 326 N 
3C 338 
3C 386 
3C 405 
3C 424 
3C 433 
3C 436 
Cen A 



0.64(0.11) 

< 4. 
0.55(0.11) 

< 0.4 
0.42(0.14) 

1.6 (0.2) 
0.12(0.06) 

< 0.12 

< 0.12 
0.30(0.06) 

< 0.2 

< 0.2 

< 1.0 
0.45(0.07) 

2.0 (0.4) 
0.48(0.16) 
10. (6.) 



1.36(0.58) 

25. (2.) 
0.41(0.22) 
0.27(0.11) 
0.55(0.29) 

5.3 (0.4) 
0.36(0.07) 

1.0 (0.1) 
0.44(0.07) 
0.69(0.06) 
0.34(0.06) 
0.43(0.06) 

1.5(0.4) 
1.12(0.06) 

1.4 (0.2) 
0.46(0.07) 

57. (9.) 



0.59(0.15) 
6.7 (0.5) 
0.23(0.09) 
0.57(0.08) 
0.27(0.15) 
1.94(0.07) 
0.08(0.04) 
0.30(0.04) 
0.25(0.04) 
0.41(0.04) 
0.32(0.10) 
0.09(0.04) 

< 1.9 
0.21(0.03) 
0.58(0.17) 
0.31(0.12) 
27. (4.) 



0.60(0.17) 

43. (5.) 
0.55(0.08) 
0.47(0.20) 
1.4 (0.2) 
3.25(0.07) 
0.48(0.04) 
0.50(0.04) 
0.44(0.05) 
1.26(0.05) 
0.81(0.12) 
0.46(0.05) 

4.1 (0.5) 
0.43(0.04) 

1.2 (0.2) 
0.28(0.04) 

9.3(3.7) 



0.76(0.27) 
8.2 (1.1) 
0.42(0.15) 

< 0.3 

< 0.6 
1.2 (0.1) 

0.19(0.06) 

< 0.11 
0.25(0.07) 
0.31(0.06) 

< 0.3 

< 0.11 

< 3. 
0.15(0.04) 

< 0.6 

< 0.10 

< 13. 



0.80(0.18) 

15. (1.) 
0.23(0.05) 
0.81(0.26) 
1.0 (0.4) 
2.8 (0.2) 

< 0.3 
0.28(0.15) 
0.58(0.11) 
0.46(0.22) 
0.20(0.10) 
0.37(0.10) 

3.5 (0.9) 
0.17(0.07) 

1.6 (0.5) 
0.22(0.10) 

12. (3.) 



< 0.9 
2.9 (0.9) 

0.37(0.19) 
1.6 (0.4) 
0.31(0.11) 

< 0.14 
0.27(0.09) 
0.44(0.11) 
0.25(0.09) 

< 0.2 
0.30(0.06) 

1.6 (0.7) 

< 0.18 

< 0.5 
0.41(0.12) 

< 8. 



< 0.4 
8.8 (0.9) 

0.16(0.05) 

< 0.7 

< 0.6 

1.3 (0.1) 
0.29(0.09) 

< 0.3 
0.17(0.09) 
0.29(0.09) 

< 0.2 
0.18(0.07) 

4.4 (0.8) 
0.28(0.05) 

1.2 (0.3) 
0.21(0.10) 
27. (4.) 



Note. — H 2 pure rotational line fluxes (10~ 14 erg s -1 cm -2 ) or 2<r upper limits measured with Spitzer IRS, 
with la uncertainties in parentheses. The 3C 218 H 2 S(6) line was unmeasurable owing to a large instrumental 
noise spike. 



TABLE 5 
Radio MOHEG PAH Features 



Source 


6.2 fim 


7.7 fim 


8.6 /urn 


11.3 fim 


12.0 /mi 


12.6 (im 


17 [im 


3C 31 


22.7(0.5) 


108.(2.) 


20.0(0.8) 


50.5(0.5) 


16.4(0.7) 


27.3(0.7) 


28.4(1.0) 


3C 84 


20. (3.) 


167.(12.) 


19. (3.) 


82. (2.) 


61. (2.) 


34. (3.) 


38. (10.) 


3C 218 


3.5(0.2) 


14.(2.) 


4.0(0.5) 


8.8(0.4) 


< 0.8 


3.9(0.5) 


2.8(1.3) 


3C 270 


5.2(0.6) 


17.(2.) 


< 0.9 


7.0(0.2) 


5.7(0.4) 


4.1(0.4) 


< 0.5 


3C 272.1 


16. (1.) 


21.(4.) 


< 6. 


28.9(0.5) 


5.2(0.6) 


11.4(0.7) 


22.1(1.4) 


3C 293 


12.0(0.6) 


50.(2.) 


21.9(0.6) 


42.7(0.5) 


14.4(0.5) 


23.6(0.6) 


26.2(1.5) 


3C 310 


0.56(0.28) 


< 1.4 


< 0.3 


0.42(0.12) 


0.29(0.13) 


0.15(0.06) 


< 0.3 


3C 315 


1.6(0.4) 


3.6(1.0) 


0.52(0.23) 


3.0(0.2) 


0.87(0.22) 


1.9(0.3) 


4.8(0.6) 


3C 317 


< 1.3 


< 2. 


< 0.8 


0.85(0.13) 


0.92(0.20) 


1.3(0.2) 


1.5(0.4) 


3C 326 N 


< 0.5 


< 0.6 


< 0.4 


0.57(0.12) 


< 0.3 


< 0.2 


< 0.7 


3C 338 


0.70(0.32) 


< 3. 


< 1.0 


1.9(0.3) 


1.0(0.4) 


< 0.8 


< 0.8 


3C 386 


0.44(0.21) 


< 2. 


< 0.6 


1.1(0.1) 


< 0.4 


0.72(0.18) 


< 0.7 


3C 405 


< 6. 


< 22. 


< 6. 


< 11. 


< 14. 


< 12. 


< 14. 


3C 424 


< 0.3 


< 0.6 


< 0.2 


0.42(0.08) 


0.57(0.11) 


0.87(0.15) 


< 0.6 


3C 433 


2.5(1.1) 


8. (4.) 


9. (1.) 


3.6(0.8) 


9.5(1.0) 


13. (1.) 


32. (2.) 


3C 436 


< 0.7 


3.2(1.6) 


< 0.8 


1.1(0.4) 


< 1.4 


2.1(0.7) 


< 1.6 


Cen A 


198.(17.) 


1052.(92.) 


459.(32.) 


1557.(26.) 


848.(24.) 


655.(27.) 


2668.(73.) 



Note. — Flux (10 14 erg s 1 cm 2 ) or 2a upper limit measured with Spitzer IRS and PAHFIT, 
with la uncertainties in parentheses. 
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TABLE 6 

H 2 Upper Limits or Single-Line Detections and PAH Features 



Source 


H 2 S(0) 


H 2 S(l) 


H 2 S(2) 


H 2 S(3) 


PAH 6.2 


PAH 7.7 fim 


PAH 11.3 fim 


3C 15 


< 0.4 


< 0.4 


< 0.10 


< 0.09 


< 0.9 


< 2. 


0.57 (0.13) 


3C 29 


0.55(0.15) 


< 0.7 


< 0.2 


< 0.4 


< 0.9 


< 3. 


< 0.44 


3C 66B 


< 0.12 


0.30(0.07) 


< 0.10 


< 0.11 


< 0.5 


< 0.7 


0.73 (0.10) 


3C 76.1 


< 0.2 


0.45(0.09) 


< 0.12 


< 0.06 


< 1.3 


< 3. 


1.48 (0.11) 


3C 83.1 


< 0.14 


< 0.14 


< 0.11 


0.25(0.06) 


1.5 (0.5) 


3.2 (1.6) 


2.7 (0.2) 


3C 120 


< 2.2 


2.74(0.75) 


< 1.5 


< 1.0 


7.9 (2.0) 


117. (14.) 


33. (2.) 


3C 129 


< 0.3 


< 0.19 


< 0.16 


0.31(0.08) 


< 0.9 


< 1.5 


2.9 (0.2) 


3C 189 


0.23(0.08) 


< 0.4 


< 0.11 


< 0.3 


< 1.3 


< 1.7 


< 0.26 


3C 264 


< 0.3 


0.28(0.12) a 


< 0.11 


0.20(0.05) 


7.6 (0.5) 


14.9 (1.7) 


6.2 (0.2) 


3C 274 


< 0.7 


< 1.1 


< 1.6 


2.92(0.62) 


< 7. 


< 28. 


< 3. 


3C 318.1 


< 0.14 


< 0.16 


< 0.06 


< 0.06 


< 0.7 


< 1.2 


0.31 (0.07) 


3C 348 


< 0.4 


< 0.5 


< 0.4 


< 0.5 


< 2. 


< 3. 


< 1.6 


3C 403.1 


< 0.2 


< 0.16 


< 0.08 


0.10(0.03) 


< 0.4 


< 1.1 


0.65 (0.10) 


3C 449 


< 0.09 


< 0.4 


< 0.08 


< 0.04 


< 0.2 


< 0.6 


0.24 (0.07) 


3C 465 


< 0.18 


< 0.2 


< 0.07 


< 0.10 


< 0.4 


< 1.1 


0.58 (0.10) 


NGC 6251 


< 0.3 


< 0.3 


< 0.18 


< 0.2 


1.6 (0.2) 


6.6 (1.0) 


2.7 (0.2) 


IC 4296 


< 0.3 


< 0.3 


< 0.2 


0.78(0.15) 


< 3. 


5.6 (1.0) 


2.6 (0.3) 


3C 17 


< 0.4 


< 0.5 


< 0.5 


< 0.2 


< 1.2 


< 2. 


1.6 (0.3) 


3C 28 


< 0.3 


< 0.4 


< 0.4 


0.22(0.04) 


< 1.8 


< 1.2 


5.8 (0.6) 


3C 33 


< 0.7 


< 0.5 


< 0.19 


0.25(0.09) 


< 2. 


< 6. 


5.2 (0.4) 


3C 61.1 


< 0.3 


0.37(0.10) 


< 0.12 


< 0.11 


< 1.4 


< 3. 


< 0.5 


3C 111 


< 1.2 


< 1.4 


< 1.4 


< 0.9 


< 6. 


< 28. 


< 4. 


3C 123 


< 0.4 


< 1.7 


< 0.7 


< 0.3 


< 1.7 


< 0.9 


< 1.3 


3C 192 


< 0.2 


< 0.3 


< 0.16 


< 0.13 


< 1.1 


< 1.7 


< 0.5 


3C 219 


< 0.3 


0.45(0.09) 


< 0.3 


< 0.16 


1.9 (0.6) 


< 3. 


< 0.5 


3C 234 


< 1.9 


b 


< 0.9 


< 0.4 


4.6 (1.6) 


9.2 (4.4) 


< 2. 


3C 319 


< 0.3 


< 0.2 


< 0.3 


< 0.12 


< 0.7 


< 1.0 


0.31 (0.12) 


3C 321 


< 1.8 


< 1.5 


< 1.1 


1.44(0.39) 


< 7. 


61. (11.) 


6. (2.) 


3C 346 


< 0.18 


< 0.3 


< 0.10 


< 0.15 


< 0.4 


< 1.1 


< 0.3 


3C 381 


< 1.5 


< 1.2 


< 0.3 


< 0.2 


< 1.0 


3.8 (1.0) 


1.2 (0.2) 


3C 382 


< 0.4 


< 0.6 


< 0.4 


< 0.7 


< 2. 


< 6. 


< 2. 


3C 388 


< 0.2 


< 0.3 


< 0.17 


< 0.2 


< 0.9 


< 1.1 


< 0.6 


3C 390.3 


< 0.9 


< 0.6 


< 0.3 


0.24(0.09) 


< 1.8 


14.1 (1.5) 


< 2. 


3C 401 


< 0.4 


< 0.2 


< 0.11 


< 0.13 


< 0.7 


< 1.1 


< 0.4 


3C 445 


< 1.8 


< 1.9 


< 1.6 


< 0.8 


< 4. 


7.5 (3.6) 


< 3. 


3C 452 


< 0.5 


< 0.3 


< 0.3 


< 0.18 


< 1.5 


3.2 (1.7) 


< 1.2 


3C 459 


< 2.6 


< 1.1 


0.94(0.34) 


< 0.4 


9.6 (2.1) 


15.0 (1.3) 


9.6 (0.6) 


Pic A 


< 0.3 


< 0.7 


< 0.7 


< 0.7 


< 6. 


< 8. 


< 1.6 



Note. — H 2 pure rotational line and PAH fluxes (10 14 erg s 1 cm 2 ) or 2cr upper limits measured with 
Spitzer IRS, with lcr uncertainties in parentheses. 

a The 3C 264 H 2 S(l) line is only a marginal detection (2.3<r), so this source does not meet our H 2 detection 
criterion of 2 lines detected at > 2.5o\ 

b The 3C 234 H 2 S(l) line was unmeasurable owing to a large instrumental noise spike. 



TABLE 7 

Radio MOHEG Forbidden Emission Lines 



Source 


[Fe II] 


[Ar n] 


[Ar in] 


[S iv] 


[Ne II] 


[CI n] 


[Nc m] 


[S in] 


[O iv] 


[S in] 


[Si n] 




5.34 (im 


6.99 fim 


8.99 Atm 


10.51 //m 


12.81 Atm 


14.37 lira 


15.56 (im 


18.71 lira 


25.89 /zm 


33.48 Atm 


34.82 fim 


3C 31 


0.48(0.18) 


1.3 (0.2) 


< 0.3 


0.51(0.18) 


3.0 (0.2) 


< 0.4 


1.5 (0.2) 


0.73(0.17) 


1.1 (0.2) 


0.84(0.15) 


2.4 (0.3) 


3C 84 


10. (1.) 


12.4(0.9) 


2.5 (0.8) 


1.4 (0.6) 


43.3 (0.5) 


6.8 (1.7) 


26. (2.) 


13. (2.) 


< 7. 


23. (2.) 


62.3 (0.6) 


3C 218 


0.14(0.05) 


0.46(0.07) 


< 0.3 


< 0.2 


1.78(0.08) 


< 0.5 


0.78(0.17) 


< 0.5 


1.0 (0.2) 


0.93(0.26) 




3C 270 


< 0.8 


2.4 (0.3) 


0.23(0.10) 


0.33(0.08) 


4.20(0.08) 


0.58(0.20) 


2.3 (0.1) 


0.87(0.11) 


1.3 (0.3) 


1.3 (0.2) 


2.9 (0.1) 


3C 272.1 


< 1.5 


3.8 (0.6) 


< 0.6 


1.0 (0.2) 


5.6 (0.2) 


0.44(0.20) 


5.7 (0.2) 


3.4 (0.2) 


1.4 (0.2) 


4.1 (0.1) 


11.7 (0.3) 


3C 293 


1.3 (0.2) 


1.9 (0.2) 


0.61(0.09) 


< 0.3 


5.24(0.08) 


< 0.3 


2.6 (0.1) 


1.2 (0.2) 


1.5 (0.1) 


2.8 (0.3) 




3C 310 


< 0.2 


< 0.3 


0.09(0.04) 


< 0.08 


0.28(0.04) 


< 0.2 


0.47(0.06) 


< 0.2 


< 0.2 


1.1 (0.2) 




3C 315 


0.28(0.11) 


0.53(0.14) 


0.12(0.04) 


0.10(0.04) 


0.29(0.05) 


< 0.2 


0.30(0.09) 


< 0.3 


0.24(0.08) 






3C 317 


< 0.3 


0.65(0.13) 


0.21(0.05) 


< 0.08 


0.88(0.05) 


0.19(0.07) 


0.43(0.07) 


0.37(0.07) 


0.26(0.08) 


< 0.3 


1.3 (0.2) 


3C 326 N 


0.23(0.08) 


< 0.3 


0.30(0.05) 


0.09(0.04) 


0.44(0.05) 


0.17(0.05) 


< 0.2 


< 0.2 


0.12(0.05) 






3C 338 


< 0.3 


< 0.4 


< 0.2 


0.22(0.11) 


0.66(0.11) 


< 0.2 


0.27(0.07) 


< 0.3 


0.34(0.10) 


0.58(0.16) 


1.2 (0.3) 


3C 386 


0.43(0.08) 


< 0.2 


< 0.1 


< 0.1 


0.39(0.04) 


0.13(0.06) 


0.61(0.05) 


0.39(0.10) 


< 0.3 


0.96(0.21) 


2.9 (0.2) 


3C 405 


2.3 (0.9) 


11.6 (0.8) 


4.4 (0.7) 


16.2 (0.7) 


21.7 (0.7) 


5.4 (0.4) 


47.9 (0.5) 


24.2 (0.6) 


78.5 (0.7) 


29. (1.) 




3C 424 


0.15(0.05) 


< 0.07 


< 0.09 


< 0.07 


0.25(0.08) 


0.12(0.06) 


0.16(0.05) 


< 0.2 


0.18(0.06) 






3C 433 


< 0.6 


1.3 (0.4) 


1.5 (0.2) 


1.7 (0.2) 


2.6 (0.2) 


< 0.9 


5.7 (0.2) 


0.89(0.20) 


7.0 (0.4) 






3C 436 


< 0.2 


< 0.1 


0.12(0.05) 


0.12(0.04) 


0.29(0.10) 


< 0.3 


0.85(0.10) 


0.28(0.11) 


0.93(0.12) 






Cen A 


< 18. 


67. (4.) 


10.9 (5.2) 


28.(4.) 


213.(4.) 


37. (12.) 


233.(18.) 


87. (9.) 


155. (18.) 


264.(11.) 


561.(14.) 



Note. — Flux (10 14 erg s 1 cm 2 ) or 2<r upper limit measured with Spitzer IRS, with la uncertainties in parentheses. 



Molecular Hydrogen in Radio Galaxies 



TABLE 8 

HlGH-IONIZATION FORBIDDEN EMISSION LINES 

Source [Ne vi] 7.65 [Ne v] 14.3 [Ne v] 24.31 

3C 405 13. (1.) 20.5(0.5) 24.3(0.6) 

3C 433 2.3(0.4) 0.8(0.2) 0.9(0.3) 

3C 436 0.23(0.05) < 0.3 < 0.2 

Cen A 29.(14.) 37.(11.) 65.(31.) 

Note. — Flux (10 -14 erg s _1 cm -2 ) measured 
with Spitzer IRS, with la uncertainties in parenthe- 
ses. 



TABLE 9 

Radio MOHEG Continuum and Line Luminosities 



Source logi/L„(178 MHz) logi/L„(24 /xm,L Q ) logL(H 2 ) a logL x (2-10 keV) b 



3C 31 


40.33 


8.63 


40.32 


40.67 


3C 84 


40.92 


10.84 


41.75 


42.91 


3C 218 


42.46 


9.34 


41.10 


41.69 


3C 270 


40.07 


8.23 


39.26 


41.08 


3C 272.1 


39.02 


7.42 


38.86 


39.34 


3C 293 


41.07 


9.74 


41.76 


42.78 


3C 310 


41.87 


8.64 


40.85 




3C 315 


42.02 


9.43 


41.76 


41.68 c 


3C 317 


41.42 


8.61 


40.53 


41.30 


3C 326 N 


41.90 


8.57 


41.73 


40.63 c 


3C 338 


41.28 


8.49 


40.54 


40.30 


3C 386 


40.48 


7.75 


39.87 


39.75 c 


3C 405 


44.10 


11.31 


41.75 


44.28 


3C 424 


42.08 


9.53 


41.97 




3C 433 


42.46 


11.15 


42.13 


43.90 c 


3C 436 


42.67 


10.32 


42.31 


43.53 


Cen A 


40.12 


9.07 


39.88 


41.70 



a Observed H2 rotational line luminosity (erg s -1 ), summed over the 0-0 S(0)-S(3) tran- 
sitions (Table 4). Corrections of 6-33% are made for undetected emission lines, based on 
H2 multi-temperature models (§5). 

b Un a bsorbed 2-10 keV nuclea r X-ray luminosit i es (er g s -1 ) from I Evans et al. I 
(200 61). iBalmaverde et al. I 120061), [ Hardcastlc et al. I i200B ), fMcNama ra et al. I 112000). 
IPellegrinil (120051) . and lUeda et al. I (120051) 

Unabsorbed 2-10 keV nuclear X-ray luminosities (erg s — 1 ) measured with Chandra 
(Table 10). 



TABLE 10 

Chandra X-ray Data and AGN Spectral Models 



Source 


t (ks) 


chip 


model 


ix a 


Notes 


3C 315 


7.7 


S3 


PL 


4.8 x 10 41 


S/N too low to detect heavily abs. emission 


3C 326 N 


27.5 


13 


PL 


4.3 x 10 40 


Low S/N 


3C 386 


29.1 


13 


PL 


5.6 x 10 39 




3C 433 


37.2 


S3 


PL + abs. PL 


8.0 x 10 43 


See Miller & Brandt (2009) 


Zw 3146 


45.6 


13 


PL 


8.2 x 10 42 


Likely contamination from thermal gas 



a Unabsorbed 2-10 keV nuclear X-ray luminosities (erg s 1 ), measured with Chandra. 



Ogle et al. 



TABLE 11 
H 2 Model Parameters 



Source 


T(K) a 


Ortho/Para b 


iV(H 2 ) c 






M(H 2 ) d 


Ltot(H 2 ,rot) c 


op qi 


100. ( 0.) 


l.Oo I 


o.O-Cj-tZI 


z 


QP_i_nfi 
oHj-|-Uo 


( 1 fiP_i_nfi q 1 p 1 n«\ 
l.orj+Uo,o. 1Hj-|-Uo ) 


QP 1 QQ 

o.o-cj-f-oy 




250. 


( 60.) 


z.yoo 


9 QT7 OH 





HP I flfi 

Ui!j-r-UD 


i~\ fiP_i_nfi 9 npj_n7\ 
( l.o.cj-|-UO,z.U-Cj-t-U ' ) 


1 fip 1 /i n 

1.0.CjH-4U 




1140. 


150.) 


^ nnn 


o ^p_i_i 7 


9 
Z 




(\' ^"P-LnA a nTr_i_nA^ 

( 1 . Oll;-|-U i t,4.Ull; - rU4 J 


9 ^pj_An 

Z.Oll; _ r4U 


OP Q/t 


150. 


( 20-) 


9 ARA 
Z.4D4 








<^orij-T"Uo 


9 npj_Ai 

Z.U-Cj _ r41 




580. 


( 30.) 


^ nnn 

O.UUU 


1 7P_L_9n 
1. 1 JZj-tZU 


A 
'1 


7P_i_nfi 


^.Ullj-l-UOjO.Olij-rUO J 


Q OP HI 
0.0-Cj-t41 




1500. 


( 0.) 


^ nnn 

O.UUU 


^ 9P-l_1 £ 




y 




?s np_j_nA 1 nTTj-n^^ 

( O.Ull;-|-U'y:, l.U-Cj-rUO ) 


9 9P_LA1 
Z.ZHj-t-41 


0\y Zlo 


100. 


( !0-) 


1.004 


£ 9P_l_91 
O.Z-Ej-^ZI 


9 
z 


nP-uno 


(~\ 9P_Lno ^ ^p_Lno\ 
( i.zi!j-|-uy,o.O-cj-ruy ) 


7 9P_LAn 
( .Zll;-t-4U 




540. 


( 70.) 


^ nnn 

O.UUU 


A 7P-L1 £ 


1 

1 


9P_i_nfi 


(7 nP-L_n^ 9 nP-Lfi^i 

. ( .UH;-|-UO,Z.U-Cj-rUO ) 


Q AP_LAn 
0.41L;-f-4U 




1500. 


(0.) 


^ nnn 

O.UUU 


^ QP_i_1 fi 


^ 




(a nP-unA 9 iP-unA^ 


a nP-i-An 


0^ z / u 


100 


(0.) 


1 ^£7 
l.Oo / 


\Orij-rZi 






^^Pj_n7 
<^0-Cj-|-U i 


D. 1 rjTOo 




400. 


( 70.) 


O OOO 

z.yyy 


Q 1 P_l_1 O 
0. l-Cj-f-iy 


1 

1 


fit? 1 nc; 
OHi-|-UO 


/ fi QP_l_n/1 9 fiP_LO^\ 
( O.ol!j-r-U4,Z.O-Cj-)-U0 _J 


9 fiP_LQO 

z.oiij-)-oy 




1500. 


(0.) 


^ nnn 

O.UUU 


Q IP 11 7 


1 

1 


7T? _un^ 


?i AP_i_n^ 9 iP-Ln'i^ 

1 .4ij-t-UO,Z. lUi-f-UO J 


o.yiij-roy 


op 979 1 
0^ Z / Z. 1 


100 


(0.) 


1 ^S7 


7P_I_01 


7 
< 


op 1 

yHj-hUD 


/ ^ QP I nfi 1 0P_L07\ 

( 0.oi!j-r-UD,l.zli;-)-U 1 ) 


1 OP I QO 

i.yn;-)-oo 




290. 


180.) 


z.yoo 


O VT7 I 1 O 


Q 
O 


OT? 1 r\A 


( O.4i1/+Uo,o.U±!j+U0 ) 


1.0-Cj+OO 




1150. 


120.) 


q nnn 
o.UUU 




9 
Z 


OT? 1 flQ 

ZHjH-Uo 


?1 c;P_l_nQ Q OPJ_nQ\ 
( 1.0l!j-hUo,O.Zli;-)-Uo ) 


9 flP _L ffl 

z.uiijH-oy 


9f 9QQ 

oLj zyo 


100 


(0.) 


1 ^87 
l.Oo 1 


9 1 P-L99 

z. iej-tzz 


O 


(^p_i_no 


?q 1 p_i_no a 9P_i_no\ 
( o. irj-|-uy ,4.zii;-t-uy ) 


1 OP-LAI 
l.UI-; _ r41 




278 


(9.) 


9 osn 
z.you 


A 8P_l_9n 

4.0-CJ^ZU 


Q 

O 


OCjtU 1 


(7 np_j--n7 q 7Pj_n7^ 


A 1 P-LA1 
4. 1Jzj-|-41 




1040. 


( 30.) 


*z nnn 

O.UUU 


o.ZHj-r lo 







?a sp_i_n^ <^ 9P_Ln r i^ 

( 4.oil;-|-UO,O.Zll;-t-UO J 


A 1 P-LA1 
4. lrj-|-41 


op Tin 


140. 


( 20-) 


9 Q/in 

Z.o4U 


7 1 t? on 
< . lUi+ZU 


1 

1 


7P_i_ns 


/o np_i_n7 q opj_na\ 
( y.Ui!j-r-U 1 ,o.Zli;-)-Uo J 


o.4ii;-|-4U 




640. 


160.) 


q nnn 
o.UUU 


T 71? _L 1 S 

1. ( Hj-f-lo 


A 
1 


i P 1 n£ 
Irj-hUO 


?o ip 1 nt; a npj_nc;\ 
^Z. 1i1/H-UO,o.U-Cj-)-UO J 


7P_l /in 

O. / HjH-4U 




1500 


( 0.) 


*z nnn 

O.UUU 


1 1 P_l_1 7 
1. lHj-p 1 ( 


9 
Z 


717 _unA 


?i op_i_nA q sP-LnA^ 
( 1 .yi!j-|-U4,o.oi!j _ rU4 / 


7 AP_l_An 

( .4H;-t-4U 


On QIC 


160. 


( 40-) 


Z.DUo 


,- /i p_i_on 
<^4ll/-t-ZU 






/QT?_Lflfi 

<.otii-\-\Jo 


T QP_L/1T 
1 .oHit^I 




330. 


( 40-) 


9 QQA 

z.yy4 


o.O-cj-t iy 


O 


1 p_i_n7 


(A AP-L_n7 7 nP 1 ^n7 , \ 

, 1 . 4lZj _ pU I;'. UJ_j 1 U 1 ) 


O.olij-|-41 




1500 


(0.) 


^ nnn 

O.UUU 


I.OJZj-T 1 i 


]^ 


ziP-un^ 


( 1.1HjA^UO,1. 1 Hi^UO ) 


A fiP-UA1 


op 01 7 


100 


(0.) 


1 ^87 

l.OO 1 


<2E+21 






/9P 1 nc 

<^ZHj-)^UO 


A ^P4_^Q 

^.O-Cj^oy 




300. 


( 40.) 


9 Q&& 
z.yoo 


Q AP-Ul Q 
o.*±Hi^ iy 


Q 
O 


^P-unfi 

OHjh^UU 


("9 9P_i_nfi ^ ^P-unR^ 

( Z.ZHjA^UU,0.011;T^UIJ j 


9 9P_i_An 




1220. 


150.) 


*z nnn 

O.UUU 


A ^Pj-1 7 
4.0H/-r 1 ' 


1 




(1 np_i_nA « 7PJ_nA^ 

( 0.Ull;-|-U4,0. 1 H;-t-U4 ^ 


^ 7Pj_An 


of Q9P, NT 
OZO IN 


110. 


( !0-) 


1 7M 
1.(01 


'i 7P_I_91 
O. 1 rij-t-Zl 


9 
Z 


9P_i_no 


?i QP_i_no ^ sp_i_no\ 
( i.oii;-|-uy,o.o-cj-ruy ^ 


1 1P-LA1 
1. lH;-t _ 41 




430. 


( 60.) 


q nnn 
o.UUU 


1 OP_i_ 1 


1 

1 


op_i_n7 


/7 7P_l_nfi 1 7P_L07\ 
( ( . ( Hj-|-UD,1. 1 Hj-)-U i j 


0.yHj-)-41 




1110. 


250.) 


*z nnn 

O.UUU 


O.DIjt 1 ' 


Q 
O 


41j-)-UO 


?i 9P_i_n^ q 7Pj_n r i^ 
1 . zi!j-|-uo,y. ( i-j i uo ) 


0.oj_j-(-'11 


of QQQ 


100 


(0.) 


1 ^£7 
l.Oo 1 


^ op 1 1 ni 






^9P_i_ns 


u.O-cj-roy 




480. 


( 60.) 


^ nnn 

O.UUU 


Q AP-Ul & 
y .'ilLjn^ 10 


7 


fiP_un^ 

UHjh^UO 


( 0.1HjA^UO,1.111/t^UU ) 


^ ^P-uAn 

O.O-Cj^^U 




1450. 


590.) 


^ nnn 

O.UUU 


, - 1 p 1 1 7 






1 p 1 HA 


9 9P_LAn 
Z.Z_Cj-t"4U 


op OOP 
OO OOU 


130. 


( 20.) 


9 1 9 A 


.r-" 9 P _U 9 1 
^ZEj-pZl 






//IP 1 H7 


o.ynj-roy 




430. 


150.) 


q nnn 
o.UUU 


^ OP I 1 Q 
O.Zil/-|- lo 


1 


at? 1 nc; 
4rj-|-U0 


( Z.oi!j-r-U4,0.4lL;-hUO ) 


Q fiP 1 QO 

o.0-cj-)-oy 




1500. 


(0.) 


q nnn 
o.UUU 


O 117x17 


O 


tr T7> 1 HQ 

orj+Uo 


( 4.Oi!jH-Uo,O.0i1i+Uo J 


1 QT7 I /in 


op /inf; 


100 


(0.) 


l.Oo / 


^ OP_LOO 
<^ ZHj-f-ZZ 






/ip 1 no 

■C4Hj-)-uy 


1 .Zll/-f-41 




460. 


( 70.) 


nnn 

O.UUU 


^t.onj-riy 


1 

1 


9P_i_n7 
zii/-|-u ( 


/"r 8P_i_rifi 9 ip_i_n7'\ 

( O.OI1;-|-UO,Z. llL;-t-U ( ) 


A RP-I-A1 
4.0I-; _ r41 




1500 


(0.) 


*z nnn 

O.UUU 


I.o-Cj-t lo 


1 


OCj-|-UO 


?^ QP-Ln^ ^ ^P-Ln^^ 

( 0.yH;-|-UO,0.0-Cj-rUO j 


1 ^P_LA9 
1.0H/-T4Z 


3C 424 


132 


(8.) 


2.237 


2.9E+21 


3 


3E+09 


(2.4E+09.4.4E+09) 


6.5E+41 




360. 


( 60.) 


2.997 


1.6E+19 


1 


8E+07 


(8.5E+06.4.0E+07) 


3.4E+41 




1340. 


230.) 


3.000 


2.0E+17 


2 


2E+05 


(1.2E+05.4.3E+05) 


5.4E+41 


3C 433 


100 


(0.) 


1.587 


3.1E+22 


2 


3E+10 


(1.9E+10,2.8E+10) 


8.3E+41 




350. 


( 80.) 


2.996 


4.2E+19 


3 


2E+07 


(1.4E+07.7.4E+07) 


4.7E+41 




1300. 


210.) 


3.000 


9.3E+17 


7 


1E+05 


(3.9E+05.1.3E+06) 


1.4E+42 


3C 436 


100 


(0.) 


1.587 


6.2E+21 


1 


6E+10 


(1.1E+10.2.4E+10) 


8.6E+41 




300. 


( 40.) 


2.987 


2.9E+19 


7 


6E+07 


(4.1E+07,1.4E+08) 


8.6E+41 




1500 


(0.) 


3.000 


2.0E+17 


5 


2E+05 


(4.3E+05.6.4E+05) 


2.5E+42 


Ccn A 


100 


(0.) 


1.587 


1.1E+23 


3 


4E+07 


(1.9E+07.5.9E+07) 


8.2E+38 




290. ( 30.) 


2.986 


5.0E+21 


1 


6E+06 


(l.lE+06,2.3E+06) 


8.0E+39 




1500. ( 0.) 


3.000 


9.6E+18 


3 


OE+03 


(2.7E+03,3.4E+03) 


6.7E+39 



a H 2 model component temperatures. A temperature range of 100 < T < 1500 K was enforced in 
fitting the models. Temperatures pegged at cither end of this range have indeterminate uncertainties. 
b Ratio of ortho (odd J) to para (even J) H 2 , fit self-consistently to match the thermal equilibrium 
value. 

c H 2 model column densities (cm -2 ), assuming a source size of 3'.'7 X 3'/7, matching the SL slit width. 
d H 2 model masses and lcr mass ranges (Mq). 

e H 2 model luminosities (erg s _1 ), summed over all pure-rotational transitions. Note that not all of 
these transitions are directly observed by Spitzer. 



Molecular Hydrogen in Radio Galaxies 



TABLE 12 
Radio MOHEG H 2 and Dust Masses 



oourcG 


lvi (Wcii hi r±2 ) 


ivi i coiu JH2 / 


widrrn / cold 


Ref d 


A//", e 
J «dust 


3C 31 


2.4E+08 


1.3E+09 


0.18 


4 


2.0E+05 


3C 84 


4.8E+06-8E+08 


2.2E+09 


< 0.36 


5 


4.0E+05 


3C 218 


2.0E+09 


2.0E+09 


1.0 


2 




3C 270 


1.6E+05-3E+07 


<9E+07 




8 


1.0E+03 


3C 272.1 


7.9E+06 


3.7E+06: 


2.1: 


8 


4.0E+03 


3C 293 


3.7E+09 


2.3E+10 


0.16 


1 


2.0E+06 


3C 310 


1.7E+08 










3C 315 


3.1E+07-3E+08 








1.3E+05 


3C 317 


3.5E+06-2E+08 










3C 326 N 


2.2E+09 


1.7E+09 


1.3 


9 




3C 338 


7.6E+05-2E+08 


<lE+09 




2 


1.3E+05 


3C 386 


1.5E+05-4E+07 


2.2E+08 


< 0.18 


8 




3C 405 


1.2E+07-4E+09 


<lE+09 




3 




3C 424 


3.3E+09 


< 6E+09 


>0.6 


6 




3C 433 


2.3E+10 


< 9E+09 


>3 


3 


5.0E+05 


3C 436 


1.6E+10 










Cen A 


3.6E+07 


<3E+08 f 


>0.12 


7 


5.0E+05 



a Total mass {Mq) of warm H2 components from Table 11. Mass ranges are given 
for galaxies that arc undetected in H2 S(0) but detected in other H2 rotational lines, 
since this line determines our estimate of the mass of the coldest and typically most 
massive warm H2 component. 

b Cold H2 mass (Mq), estimated from CO flux density and corrected to the same 
standard Galactic CO conversion factor of 4.6 Mq/(K km s _1 pc 2 ). Upper limits 
may be underestimated in the case of broad CO lines. Tentative (low S/N) detections 
are indicated by a colon. 

c Ratio of cols. 2 and 3. 

d Cold H 2 Refs: l)IEvans et al~l (fl999l'l 2)ISalome & Combes I ||2003T) 3) lEvans et al. I 
(2005) 4) Okuda ct al. (2005) 5) Salome et al. ( 2006) 6) Saripalli & Mack (2007), 
7) IIsrael et al. I A19901) 8) IQcana Flacquer et al. I d2010t) 9) INesvadba et al. I (120101). 
e Ma ss (Mq) of clumpy dust estimated from optical absorption I IDe Koff et al. I 
2000). This could be much less than the amount of dust associated with the warm 
H2 or CO, which might be revealed with sensitive far-IR continuum measurements. 
f Cen A cold H2 mass detected within 40" aperture. 
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TABLE 13 
Star Formation Rates 



Source 


SFR(PAH 7.7) 


SFR(PAH 11.3) 


3C 31 


0.43 (0.01) 


0.78 (0.01) 


3C 84 


71 (Cl 051 


1 35 l() 031 


3C 218 


62 (0 Oil 


1 49 <() 071 


3C 270 


0.013 (0.002) 


0.021 (0.001) 


3C 272.1 


0036 (0 00071 


019 (0 0011 


3C 293 


1.47 (0.06) 


4.81 (0.06) 


3C 310 


< 0.06 


0.068 (0.020) 


3C 315 


0.67 (0.19) 


2.1 (0.1) 


3C 317 


< 0.03 


0.055 (0.009) 


3C 326 N 


< 0.07 


0.27 (0.06) 


3C 338 


< 0.04 


0.10 (0.02) 


3C 386 


< 0.008 


0.017 (0.002) 


3C 405 


< 1. 


< 1.9 


3C 424 


< 0.16 


0.42 (0.08) 


3C 433 


1.3 (0.6) 


2.2 (0.5) 


3C 436 


2.6 (1.3) 


3.5 (1.3) 


Cen A 


0.048 (0.004) 


0.27 (0.01) 



Note. — Star formation rates (Mq yr 1 ) 
are estimated from SFR(PAH 7.7)= 2.4 X 
10- 9 L(PAH7.7)/L Q and SFR(PAH 11.3)= 
9.2 X 1O- 9 L(PAH11.3)/L0, using PAH lu- 
minosities inside the Spitzer IRS SL slit. 
SFR(PAH 11.3) should be considered an up- 
per limit for sources with high (> 0.26) 11.3/7.7 
PAH ratio. 



TABLE 14 
Jet Cavity Power 



Source Cluster P(jet cavity) a L(H2)/P(jet cavity) b 



3C 84 Perseus A 1.4E+44 2.7E-03 

3C 218 Hydra A 4.3E+44 3.4E-04 

3C 272.1 M84, Virgo 1.0E+42 1.1E-03 

3C 317 A2052 1.5E+44 2.3E-04 

3C 338 A2199 2.7E+44 1.2E-04 

3C 405 Cyg A 1.3E+45 3.8E-04 



a Jet cavity power estimated from the 4pV enthalpy required 
to inflate the cavity, divided by the cavity buoyancy timescale 
URaffertv et al. 112003) . 

b Ratio of H2 pure rotational line luminosity (summed over 0-0 
S(0)-S(3)) to jet cavity power. 



